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Abstract: We construct a class of bottom-up holographic models with physics comparable 

Nf 

to the one expected from QCD in the Veneziano limit of large Nf and Nc with fixed x = j^. 
The models capture the holographic dynamics of the dilaton (dual to the YM coupling) 
and a tachyon (dual to the chiral condensate), and are parametrized by the real parameter 
X, which can take values within the range ^ < x < ^■. We analyze the saddle point 
solutions, and draw the phase diagram at zero temperature and density. The backreaction 
of flavor on the glue is fully included. We find the conformal window for x > Xc, and the 
QCD-like phase with chiral symmetry breaking at x < Xc, where the critical value Xc lies 
close to four. We also find Miransky scaling well as Efimov-like saddle points. 

By calculating the holographic /3-functions, we demonstrate the "walking" behavior of the 
coupling in the region near and below Xc- 
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1. Introduction 

Holographic techniques have been used recently in order to understand the strong dynamics 
of gauge theories. QCD is one of the obvious targets of such a program, but strongly 
coupled gauge theories may also emerge in other contexts, namely in the physics beyond 
the Standard model (non-perturbative electroweak symmetry breaking is an example), as 
well as in condensed matter contexts. 
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An interesting mild generalization of QCD, involves SU{Nc) YM coupled to Nf Dirac 
fermions transforming in the fundamental representation (that we will still call quarks). 
This is a theory, that can be studied in the {Nc,Nf) plane. 

As usual simplifications arise in the large- A'^c limit. The standard 't Hooft large- Ac 
limit, |||], lets Ac — )■ oo keeping Nj and A = Qy^Nc finite. In this limit, the effect of 
the quarks are suppressed by powers of ^ — )• 0, and therefore it corresponds to the 
"quenched" limit. In particular interesting dynamical effects as the conformal window, 
and exotic phases at finite density, driven by the presence of the quarks are not expected 
to be visible in the 't Hooft large- Ac limit. 

In ||2| Veneziano introduced an alternative large- Ac limit in which 

Ac — ^ oo , Aj oo , —J- = X fixed , A = Qym^c fixed (1.1) 

in order to make the chiral U(l) anomaly visible to leading order in the 1/Ac expansion. 
This is the large- Ac limit we will study in this paper. 

There are several interesting issues that are accessible in the Veneziano limit. 

• The "conformal window" with an IR fixed point. The window extends from x = ^- 
to smaller values of x, and includes the Banks-Zaks (BZ) weakly-coupled region as 

i. 

• The phase transition at a critical x = Xc from the conformal window to theories with 
chiral symmetry breaking in the IR. 

• A transition region near and below Xc, where the theory is expected to exhibit "walk- 
ing behavior". The theory flows towards the IR fixed point but misses it ending up 
with chiral symmetry breaking, so that the coupling constant varies slowly over a 
long range of energies. 

• New phenomena at finite density, involving color superconductivity |Q and flavor- 
color locking 

The point the point of a quantum phase transition where the theory passes 

from an IR Conformal theory (above Xc) to a theory with a non-trivial chiral condensate 
(below Xc). Such transitions were termed conformal phase transitions in and have been 
recently argued to be due to the fusion of a UV and IR fixed points Q . The scaling of the 
condensate is similar to that of the 2D Berezinskii-Kosterlitz-Thouless (BKT) transition 
Q, and is also known as Miransky scaling, |p. Several such quantum phase transitions 
have been described recently in holographic theories in [10, 11, 12|. 



The location of the lower edge of the conformal window is determined by nonpertur- 
bative dynamics. Several estimates of the value for Xc, and for the critical value of A^^ at 



finite Ac have been put forward by using different methods |13, 14, |T^, and the boundary 
of the conformal window is also being studied actively on the lattice (see, e.g., |16|). 

The walking region near Xc as well as the value of Xc, have been of interest for a while, 
due to their potential relevance for the realization of walking technicolor, [^]. Technicolor 
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has been used as a generic name for non-perturbative electroweak symmetry breaking 
mimicking the one induced by QCD, ||l^. Although non-perturbative effects in a new 
strongly coupled gauge theory can induce electroweak symmetry breaking, generating the 
quark and lepton masses is an extra problem. A new interaction (extended technicolor) is 
usually invoked to generate the requisite couplings. However their magnitude (and therefore 
the Standard model masses) are controlled by the dimension of the scalar operator that 
breaks the electroweak symmetry. In free field theory its dimension is 3, as it is a fermion 
bilinear. However, for the SM masses to have realistic values the dimension must be 
reduced to at least around 2, i.e., the anomalous dimension of the operator should be at 



least one. This is indeed expected to be generated by a "walking" theory Notice 
that coupling between the technicolor sector and the standard model may decrease the 
dimension substantially pO[|. 



There are issues that so far have made such non-perturbative approaches to be in 



apparent confiict with data, like the value of the S parameter, [21|. It was argued that 
in cases where the strongly coupled theory is near a conformal transition the S-parameter 
can be quite different potentially evading the experimental constraints^, p^]. 

There have been several bottom- up models of technicolor [24|, mostly inspired from 
the hard wall models for mesons, |2^]. Lately there have also been top-down holographic 



models of walking behavior, ||2g, 28, 12 1. They use several contexts like favored MN 



solutions, Ig^], Df — Dj pairs, [g8[ or — Dj systems, [|12 |. 

A theory that can be compared is = 1 supersymmetric QCD with Nf flavors. The 
ground states of this theory have been found by Seiberg, |2^], and we understand several 
issues associated to low energy dynamics, including the Seiberg duality. Such a theory gives 
already several important hints on the structure expected in non-supersymmetric QCD, 



pO]. Defining again x = we have the following regimes: 



At X = the theory has confinement, a mass gap and N^. distinct vacua associated 
with a spontaneous braking of the leftover R symmetry Znc- 

At < X < 1, the theory has a runaway ground state. 

At X = 1, the theory has a quantum moduli space with no singularity. This reflects 
confinement with chiral symmetry breaking. 

At X = 1 + the moduli space is classical (and singular). The theory confines, but 
there is no chiral symmetry breaking. 

Atl + -^<x<| the theory is in the non-abelian magnetic IR-free phase, with the 
magnetic gauge group SU {Nj — Nc) IR free. 

At I < X < 3, the theory flows to a CFT in the IR. Near x = 3 this is the Banks-Zaks 
region where the original theory has an IR flxed point at weak coupling. Moving to 
lower values, the coupling of the IR SU (Nc) gauge theory grows. However near x = | 



^Recent studies that extrapolate from the Banks-Zaks region suggest that the modification of the S- 
parameter near the conformal window may be modest, [E2|. 
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the dual magnetic SU{Nf — N^) is in its Banks-Zaks region, and provides a weakly 
coupled description of the IR fixed point theory. 



• At X > 3, the theory is IR free. 

The region | < x < 3 is comparable to the conformal window that is expected in (non- 
supersymmetric) QCD. Indeed, the IR coupling of the original gauge theory is becoming 
stronger as x decreases, but for x above 3/2, a new set of IR states becomes weakly 
coupled, namely the magnetic gluons and quarks. These states have been interpreted as 



the p mesons and their supersymmetric avatars, |31]. They are massless and weakly coupled 



in this region. The regime 1 + < x < ^ does not seem to have an analogue in QCD. 
The IR theory is again an IR non-abelian gauge theory, therefore trivially scale invariant 
in the far IR, but also free. 

It was suggested already in |3^, ^ that the presence of a conformal window in = 1 
supersymmetric QCD is associated with the violation of the Breitenlohner-Freedman (BF) 
bound. A recent attempt to describe related physics was done in [p4l. 



1.1 Bottom-up models for QCD in the quenched approximation 

To construct a bottom-up holographic model for QCD in the Veneziano limit^ we need to 
first understand pure YM. The simplest bottom model for pure YM in four dimensions 
is the hard- wall model, first introduced in p^]. Despite its simplicity it could capture a 
few qualitative features of the strong interaction. More sophisticated models accounted 
for the running of the YM coupling constant, incorporating therefore the dilaton into the 
gravitational action, By simply adjusting a dilaton potential they could exhibit 

many of the properties of large- A'c YM including confinement, a mas gap, asymptotic linear 
trajectories and realistic glueball spectra at zero temperature. Moreover they fared rather 
well at finite temperature, [^], and after the tuning of two phenomenological parameters 
in the dilaton potential'^, ||3^, they could agree with lattice data both at zero and finite 



temperature, |40|. The properties of IHQCD at finite temperature were further explored in 
1 41]. Alternative Einstein-dilaton models exhibiting a cross-over rather than a first order 
deconfining transition, and matching YM finite temperature dynamics were also developed 



in |42]. Such bottom up models, were used to compute transport properties of YM, like the 
bulk viscosity and the diffusion properties of heavy quarks, Backgrounds having an 
IR fixed point or a "walking" region, where the system flows close to an fixed point, were 
studied within in the IHQCD model in 45, 46 1. In these studies the fixed point was 

introduced via the input beta function, without proper modeling of the dynamics of the 
quarks, even though Nf/Nc was large. 

To go beyond YM a new ingredient is needed, namely the flavor branes. An important 
field in this context is the order parameter for chiral symmetry breaking, dual to a complex 
bifundamental field T. In the hard wall ||2^ and soft wall |47] models for mesons, such 



^We will call such a model V-QCD from now on. 

^The bottom-up Einstein-dilaton model for large-A^c YM was termed Improved Holographic QCD 
(IHQCD). 
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a field was added using a quadratic action, and chiral symmetry breaking proceeded by 
giving to such a field a vev by hand. 

In [^] it was remarked that as the flavor sector of gauge theories in string theory arises 
from D-brane-antibrane pairs, the bifundamental field T could be naturally be identified 
with the brane-antibrane tachyon field that had been studied profusely (around flat space) 
in string theory by Sen and others, |^]. The non-linear action proposed by Sen, could 
be therefore used as a well-motivated starting point in order to study the holographic 
dynamics of chiral symmetry breaking. Several general features of this approach were 



explored in [l48| , 



Chiral symmetry breaking is dynamical and is induced/controlled by the tachyon 
Dirac-Born-Infeld (DBI) action. 

Confining asymptotics of the geometry were shown to trigger chiral symmetry break- 
ing. 

A Gell-Mann-Oakes-Renner relation is generically satisfied. 

The Sen DBI tachyon action induces linear Regge trajectories or mesons. 



• The Wess-Zumino (WZ) terms of the tachyon action, computed in string theory |51]- 
[^], produce the appropriate flavor anomalies, include the axial U{1) anomaly and 
r/'-mixing, and implement a holographic version of the Coleman- Witten theorem. 

In the context above, the analysis was done in the quenched approximation: the flavor 
sector does not backreact on the metric and dilaton. Similar results were also obtained by 



considering tachyon condensation in the Sakai-Sugimoto model [54|. 



In 1 55 1 an implementation of these ideas was performed by choosing a concrete confining 



background, that is simple and asymptotically AdS. This was the Kuperstein-Sonnenschein 



background, |56], with a constant dilaton and an AdSe soliton. In this background the 



tachyon DBI action was analyzed with the following results 

• The model incorporates confinement in the sense that the quark-antiquark potential 
computed with the usual AdS/CFT prescription confines. Moreover, magnetic quarks 
are screened. 

• The string theory nature of the bulk fields dual to the quark bilinear currents is 
readily identified: they are low-lying modes living in a brane-antibrane pair. 

• Chiral symmetry breaking is realized dynamically and consistently, because of the 
tachyon dynamics. The dynamics determines the chiral condensate uniquely a s 
function of the bare quark mass. 

• The mass of the p-meson grows with increasing quark mass, or, more physically, with 
increasing pion mass. 

• By adjusting the same parameters as in QCD (Aqcd, "^ud) a good fit can be obtained 
of the light meson masses. 
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1.2 Holographic models in the Veneziano limit 

To construct V-QCD we will put together the experience from IHQCD and the tachyon 
implementation in the quenched approximation. Putting the two together we will see 
that, under reasonable assumptions, we obtain a phase diagram which is qualitatively 
in agreement with what to expect from QCD in the Veneziano limit. Moreover we will 
verify that changes of the bulk tachyon and dilaton potentials that are mild give the same 
qualitative physics. In this sense we can state confidence in our results. 
The bulk action we will consider is 

^2 



S = Sg + Sf , Sg = M^N^ J d^x^ R-^^-^ + Vg{X) 



(1.2) 



with A the 't Hooft coupling (exponential of the dilaton (f)) and the flavor action is 

Sf = -xM^N^ J d^'x Vf{X,T)^det{g^, + h{\)d^Td,T^) (1.3) 

To find the vacuum (saddle point) solution we must set the gauge fields A]^'^ to zero, as 
they are not expected to have vacuum expectation values at zero density. We also take the 
tachyon field T to be diagonal and suppressed the WZ terms as they also do not contribute 
to the vacuum solution. 

The pure glue potential Vg has been determined from previous studies, and we 
will use the same here. The tachyon potential Vj(A,T) must satisfy some basic properties, 
that are determined by the dual theory or general properties of tachyons in string theory: 
(a) To provide the proper dimension for the dual operator near the boundary (b) To 
exponentially vanish like logV/ ~ — + • • • for T — )• oo. The function h{\) captures the 
transformation from the string frame to the Einstein frame in five dimensions and will be 
chosen appropriately. 

As with IHQCD, we will arrange that the theory is logarithmically asymptotically 
AdS, and will implement the two- loop /3-function plus one-loop anomalous dimension for 
the chiral condensate. Although the geometrical picture is not expected to be reliable near 
the boundary, the renormalization group (RG) flows that emerge are reliable at least in 
the IR. The UV boundary conditions we choose can be thought of as a convenient way of 
anchoring the theory in UV. We can always define a finite cutoff and evolve the theory 
from there in the IR. 

We first analyze the fixed points of the bulk theory. Choosing a potential that imple- 
ments the Banks-Zaks fixed point, its presence exists for a range of the parameter x. We 
will make choices where this is the whole range: < x < In such a fixed point the 
dilaton is constant and the tachyon vanishes identically. We have also checked that choices 
of potential for which the fixed point exists for < x < ^ , have qualitatively similar 
physics. 

We define appropriate /3- functions for the YM coupling and the quark mass. We then 



rewrite the equations following |3€] as first order equations that specify the flow of the 
couplings, as well as non-linear first order equations that determine the /3-functions in 
terms of the potentials that appear in the bulk action ( |1.2D , (|1.3| ). 
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We can calculate the dimension of the chiral condensate in the IR fixed point theory 
from the bulk equations. We find that it decreases monotonically with x for reasonably 
chosen potentials. It crosses the value 2 at x = Xc where Xc corresponds to the end of the 
conformal window as argued in . We make the following observations which are relevant 
for technicolor studies: 

• The lower edge of the conformal window Xc lies in the vicinity of 4. Requiring the 
holographic /3-functions to match with QCD in the UV, we find that quite in general 

3.7 < xe <4.2 (1.4) 

which is in good agreement with other estimates jjl3|, 0, [l5| . 

• The fact that the dimension of the chiral condensate at the IR fixed point approaches 
two (and the anomalous dimension approaches unity) as x — >• Xc is in line with the 
standard expectation from field theory approaches [|l3| , |l9|. It is also to a large extent 
independent of the details of the model. 

It is important to stress that in the full analysis of this paper, the backreaction of the 
flavor sector on the glue sector is fully included. This is very important for the "walking" 
region, in the vicinity of x = 4, where we expect the backreaction to be important. Indeed 
we do not expect to see a Conformal Phase Transition in the quenched limit of QCD. 

Apart from x, there is a single parameter in the theory, namely where m is the UV 
value of the (common) quark mass. For each value of x, we solve the bulk equations with 
fixed sources corresponding to fixed m, Aqcd, and determine the vevs so that the solution 
is "regular" in the IR. The notion of regularity is tricky even in the case of IHQCD (pure 
glue), as there is a naked singularity in the far IR. For the dilaton this has been resolved 
in For the tachyon the notion of regularity is different and has been studied in 

detail in |5|]. 

Implementing the regularity condition in the IR and solving the equations from the 
IR to the UV (this has been done mostly numerically), there is a single parameter that 
determines the solutions as well as the UV coupling constants and vevs, and this is the 
a real number Tq controlling the value of the Tachyon in the IR. This reflects the single 
dimensionless parameter of the theory. 

For different values of x and m we find the following qualitatively different regions: 

• When Xc < X < 11/2 and m = 0, the theory fiows to an IR fixed point. The IR 
CFT is weakly coupled near x = ^ and strongly coupled in the vicinity of Xc- Chiral 
symmetry is unbroken in this regime (this is known as the conformal window). 

• When Xc < X < 11/2 and m 7^ 0, the tachyon has a non-trivial profile, and there is 
a single solution with the given source, which is "regular" in the IR. 

• When < X < Xc and m = 0, there is an infinite number of regular solutions 
with non-trivial tachyon profile, and a special solution with an identically vanishing 
tachyon and an IR fixed point. 
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• When < X < Xc and m 7^ 0, the theory has vacua with nontrivial profile for the 
tachyon. For every non-zero m, there is a finite number of regular solutions that 
grows as m approaches zero. 

In the region x < Xc where several solutions exist, there is a interesting relation between 
the IR value Tq controlling the regular solutions, and the UV parameters, namely m. This 
is determined numerically, and a relevant plot describing the relation between m and Tq 
at fixed X is in Fig. || (left). As m and —m are related by a chiral rotation by vr, we can 
take m > 0. 

The solutions are characterized by the number of times n the tachyon field changes 
sign as it evolves from the UV to the IR. For all values of m there is a single solution 
with no tachyon zeroes. In addition, for each positive n there are two solutions which exist 
within a finite range < m < where the limiting value m„ decreases with increasing 
n, and one solution for m = 0. In particular, for large enough fixed m, we find that only 
the solution without tachyon zeroes exists. 

For m 7^ 0, out of all regular solutions, the "first" one without tachyon zeroes has the 
smallest free energy. The same is true for m = 0, namely the solution with non-trivial 
tachyon without zeroes is energetically favored over the solutions with positive n as well 
as over the special solution with identically vanishing tachyon, which appears only for 
m = and would leave chiral symmetry unbroken. Therefore, chiral symmetry is broken 
for X < Xc- 

The multiplicity of regular solutions is closely related to the regime where the IR 
dimension of the chiral condensate is smaller than 2, and the associated Efimov vacua. 
They seem to be associated with the fixed point theory that here exists for all values of 
X but is not reachable by fiowing from the UV of QCD for x < Xc- On the other hand, 
the presence of a fixed-point theory in the landscape of possible theories does not seem 
necessary for the appearance of multiple saddle points. Indeed, in [55| which employed the 



quenched approximation and where no such fixed points exist, a second saddle point was 
found that provided a regular tachyon solution. It was verified however that this second 
saddle point was perturbatively unstable as meson fiuctuations were tachyonic. 

In the region just below Xc we find Miransky or BKT scaling for the chiral condensate. 
As X ^ Xc, we obtain 

^~4cDexp|^--^^=^ . (1.5) 

For X > Xc, let miR(x) be the mass of the tachyon at the IR fixed point and ^ir(x) the IR 
AdS radius. The coefficient K is then fixed as 

K= 1^ I (^-^^ 
S ["^ir^ir] x=x^ 

The construction of the holographic V-QCD model opens the road for addressing 
several interesting questions. 
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1. The calculation of the spectrum of mesons and glueballs. This is in principle a 
straightforward albeit tedious exercise, |57]. In the Veneziano limit, mixing is ex- 
pected between glueballs and mesons to leading order in 1/Nc- This will affect the 
glueball that will mix with the 0^"*" flavor-singlet o"- mesons. On the other hand 
the 2~^~^ glueballs, the 1~~ and vector mesons and the 0"' mesons do not mix, 
with the exception of the flavor singlet 0"' meson (analogous to rj') that will mix 
with the 0"' glueball due to the axial anomaly. A particularly interesting question 
here is the behavior of the mass of the lightest 0"^"*" state (the technidilaton, ||5^) as 

X — )• Xr- 



2. The structure and phase diagram of the theory at finite temperature, |59|. In the 
quenched approximation [^] the restoration of chiral symmetry was seen above the 
(first order) deconfinement transition. The expected structure is not clear here and 
several options exist. 

3. The calculation of the energy loss of heavy quarks in a quark-gluon plasma with 
non-negligible percentage of quarks. 

4. The construction of the baryon states in this theory and the calculation of their 
properties. 

5. The structure of the phase diagram at finite density and the search for exotic phases 
namely color superconductivity and color-flavor locking. 

It is plausible that the setup may provide a model for high-Tc superconductors by 
interpreting the x parameter as a "doping" parameter. The reason is that x controls the 
IR dimension of the "Cooper pair" associated with a quark-antiquark boundstate, charged 
under the axial charge. As x decreases, the IR dimension of this operator decreases, and 
the bound state becomes more and more deeply bound. At x = Xc, there is an onset of 
"axial" superconductivity (at zero temperature), that persists down to x = 0. 

At finite temperature, this picture suggests that the system might resemble the over- 
doped regime of strange metals, with x = Xc the start of the superconduction dome and 
X = the optimal doping. The connection between the value of x and doping in real 
systems may not be so far fetched as the changes in the system associated with the change 
of carriers, is accompanied by a change in the effective number of fiavors of strongly inter- 
acting effective degrees of freedom. 

The structure of this paper is as follows. In Sec. ^ we give a brief review on QCD in the 
Veneziano limit and its phase structure. In Sec. ^ we review the IHQCD model, and discuss 
the earlier results on mass spectra and the phase structure at finite temperature. Adding 
the fiavor branes is discussed in detail in Sec. ^ The V-QCD model is finally introduced 
in Sec. ^. Analysis of the model is started by studying the fixed points in Sec. ^ We 
go on transforming the equations of motion (EoMs) to equations for the holographic beta 
functions, and discuss their UV/IR asymptotics and solutions in Sec. |7[ In Sec. | we 
analyze the background, in particular how the UV expansions map to perturbation theory 
of QCD, and where the edge of the conformal window appears. We also construct and 
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present the numerical solutions for the background in the physically interesting regions. In 
Sec. we find the vacua with lowest free energy and check that they support the expected 



phase diagram. In Sec. IC we study the system near but below the conformal window, and 
show that the chiral condensate, as well as many other observables, obey the BKT scaling 
law. In particular, we check the scaling by comparing numerical results to formulas, that 



are derived analytically. Finally, we conclude and summarize the main results in Sec. |TI 
Technical details are presented in Appendices 

2. QCD in the Veneziano limit 

The conventional large- A'c limit of QCD involves a large number of colors A'c — )• oo, but 
a fixed number of fiavors, Nf — )• finite, [0. In this limit, fermion loops are suppressed, 
and the dominant diagrams are classified by the genus of the associated Riemann surface. 
Fundamentals (quarks) are associated with open strings and boundaries, and the number 
of flavors is measuring the Chan-Paton factors of the open strings. 
There is an alternative large-A'^c in QCD, in which 

Nf 

iVc-^oo , iV/->oo , -T^ = x fixed (2.1) 

This was first introduced by Veneziano in |||] in order to have the QCD axial anomaly, of 
order 0{NcNf) appear in the leading order in the large- A^c expansion. 

This alternative large- AT^ limit is very interesting in order to preserve important effects 
due to quarks. In the conventional 't Hooft limit such effects are subleading, and this is 
known as the quenched limit for flavor. Many efforts have been made in the last few years 
to consider unquenched flavor, in order to estimate the contribution of quarks to the physics 
of the quark-gluon plasma. Such efforts are summarized in the recent review, pO[| . 

The following effects are not easily visible in the conventional 't Hooft limit: 

• The "conformal window" with a non-trivial fixed point, that extends from x = ^ to 
smaller values of x. The region x — >• ^ has an IR fixed point while the theory is still 
weakly coupled, as was analyzed by Banks and Zaks, Q. 

• It is expected that at critical Xc, the conformal window will end, and for x < Xc, the 
theory will exhibit chiral symmetry breaking in the IR. This behavior is expected to 
persist down to x = 0. Above x > Xc the IR theory is CFT, at strong coupling that 
progressively becomes weak as x — )• ^. 

• Near and below Xc, there is the transition region to conventional QCD IR behavior. 
In this region the theory is expected to be "walking" , so that the theory flows towards 
the IR fixed point but misses it ending up with chiral symmetry breaking. But the 
approach to the fixed point involves a slow variation of the YM coupling constant for 
a long range of energies. This has been correlated with a nontrivial dimension for 
the quark mass operator near two, rather than three (the free field value). 

• The existence of this "walking" region makes the theory extremely interesting for 
applications to strong-couplings solutions to the hierarchy problem (technicolor). 



- 11 - 



• New phenomena are expected to appear at finite density driven by strong couphng 
and the presence of quarks. These involve color superconductivity and flavor-color 
locking [^. 

To discuss the structure expected as a function of the finite ratio x, defined in ( p.l] ) 
we write the two-loop QCD /3-function. With Nj (non-chiral) flavors in the fundamental, 
the /3-function reads 



(4^) 



11 



-Nr- 



(47r) 



34 



-N 



13 



-N 



Using the 't Hooft coupling, and setting 



Nc 



A = g'N, 



A 



X we obtain 
v3 , 



-6oA^ + 6iA^ + 0(A" 



with 



^0 = 77 



2 (11 
3 



2x) hi _ 3 (34 - 132;) 



The Banks-Zaks region is x 
point of the /3-function at 



(47r)2 
= 11/2 



?)2 
"0 



2 (11 - 2x)2 
e with e ^ 1 and positive. 



A^< 



75 



+ ••• (2.2) 



(2.3) 
(2.4) 

We obtain a fixed 
(2.5) 



which is trustable in perturbation theory, as A* can be made arbitrarily small. 

The infrared fixed point has properties that are computable in perturbation theory. 
In particular the low- lying operators consist of the conserved stress tensor, Tr\F'^\ that is 
now slightly irrelevant, and the L, R currents that are still conserved with the exception of 
the U{1) axial current that its conservation is broken by the anomaly. 

The mass operator, tpLipR now dimension slightly smaller than three, as attested 
by its perturbative anomalous dimension 

d\ogm _ oq 2 , '^1 



3 N}. 



1 



d\og fi 
1 



7 



47r^ (47r)2 



9' + 



At large Nc this becomes 

7 



ai 



2(47r) 



1) 



10 m 



2iV2 



Nc 



\+'-^iN^ 



(203 - lOx) 



X^ + 0{X^,N- 



(2.6) 
(2.7) 

(2.8) 



(47r)2'" ' 12(47r)4 

One can still perturb this theory by the [/(A''c)-invariant mass operator (assuming 
all quarks have the same mass), and the theory is expected now to flow to the trivial 
(QCD-hke) theory in the IR. 

It is believed that there is also a value Xc with < Xc < so that for x < Xc the 
theory flows to a trivial theory (with a mass gap) in the IR, with chiral symmetry breaking 
and physics isomorphic to that of standard YM. For ^ > x > Xc, the theory is expected to 
flow to a non-trivial IR fixed point, and chiral symmetry to remain unbroken, as happens in 
the BZ region. Generically, the IR theory is strongly coupled except in the region x — )■ ^ 
where the fixed point theory is weakly coupled (Banks-Zaks fixed point). For x > ^ the 
theory is IR free. 
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3. A step back: Bottom-up models for large-Nc YM 



The holographic dual of large Nc Yang-Mills theory, proposed in , is based on a five- 
dimensional Einstein-dilaton model, with the action^: 



f dl'x^ R-Ud^f + V{^) +2MInI f d^xVhK. (3.1 



Here, Mp is the five-dimensional Planck scale and Nc is the number of colors. The last 
term is the Gibbons-Hawking term, with K being the extrinsic curvature of the boundary. 
The effective five-dimensional Newton constant is G5 = l/(167rAfp A'^^), and it is small in 
the large- A^c limit. 

Of the 5D coordinates {xj, r}j=o...3) Xi are identified with the 4D space-time coor- 
dinates, whereas the radial coordinate r roughly corresponds to the 4D RG scale. We 
identify A = e* with the running 't Hooft coupling At = NcgyM ' ^° priori unknown 
multiplicative factor^, A = kXi- 

The dynamics is encoded in the dilaton potential^, V^(A). The small-A and large- 
A asymptotics of V{\) determine the solution in the UV and the IR of the geometry 
respectively. For a detailed but concise description of the UV and IR properties of the 



solutions the reader is referred to Section 2 of |38]. Here we will only mention the most 
relevant information: 

1. For small A, V'(A) is required to have a power-law expansion of the form: 

V{X) ^^{l + viX + V2X^ + ...), A^O. (3.2) 

The value at A = is constrained to be finite and positive, and sets the UV AdS scale 
i. The coefficients of the other terms in the expansion fix the /3-function coefficients 
for the running coupling X{E). If we identify the energy scale with the metric scale 
factor in the Einstein frame, denoted by below, we obtain |3(:]: 



dX , .9 , >^ , 9 , 9 207 , 

"W'SbiE"-'"' ■ '"'S"- = -4"^ + 256"' ■ 



2. For large A, confinement and the absence of bad singularities^ require: 

y(A)~A^«(logA)- A ^00, I 2/3 < Q < 2V2/3, P arbitrary 

^ ^ ^ ^ ^ ' I g = 2/3, P > 



(3.4) 



^Similar models of Einstein-dilaton gravity were proposed independently in to describe the finite 
temperature physics of large Nc YM. They differ in the UV as the dilaton corresponds to a relevant operator 
instead of the marginal case we study here. The gauge coupling e* also asymptotes to a constant instead 
of zero in such models. 

^This relation is well motivated in the UV, although it may be modified at strong coupling (see ||3^ ). 
The quantities we will calculate do not depend on the explicit relation between A and At. 

^With a slight abuse of notation we will denote V{X) the function V{^) expressed as a function of 
A = e . 

'^For a description of the notion of "bad versus good singularities" and their resolution the reader is 
referred to ]6l| ]. 
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In particular, the values Q = 2/3, P = 1/2 reproduce an asymptotically- linear glue- 
ball spectrum, ~ n, besides confinement. We will restrict ourselves to this case 
in what follows. 

In [^], the single phenomenological parameters of the potential was fixed by looking 
at the zero-temperature spectrum, i.e. by computing various glueball mass ratios and com- 
paring them to the corresponding lattice results. The masses are computed by deriving the 
effective action for the quadratic fluctuations around the background, |^2| and subsequently 
reducing the dynamics to four dimensions. 

The glueball spectrum is obtained holographically as the spectrum of normalizable 
fluctuations around the zero-temperature background. In IHQCD the relevant fields are 
the 5D metric, one scalar field (the dilaton), and one pseudoscalar field (the axion that 
is subleading in N^). As a consequence, the only normalizable fluctuations above the 
vacuum correspond to spin and spin 2 glueballs (more precisely, states with J^*^ = 
0^^,0 ^,2^+), each species containing an infinite discrete tower of excited states. 

We only compare the mass spectrum obtained in our model to the lattice results for 
the lowest 0^''',0 ^,2"*"^ glueballs and their available excited states. These are limited to 



one for each spin species, and none for the spin 2, in the study of |63], which is the one 
we use for our comparison. This provides two mass ratios in the CP-even sector and two 
in the CP-odd sector. 

The glueball masses are computed by first solving numerically Einstein's equations, 
and using the resulting metric and dilaton to setup an analogous Schrodinger problem for 
the fluctuations, |^^. The results for the parity-conserving sector are shown in Table |l|, 
and are in good agreement with lattice data for = 3. 



Table 1: Glueball Masses 





IHQCD 


iVc = 3 


Nc = oo 


mo*++/mo++ 


1.61 


1.56(11) 


1.90(17) 


m2++ /mQ++ 


1.36 


1.40(4) 


1.46(11) 



Unlike the various mass ratios, the value of any given mass in AdS-length units (e.g. 
mo++^) does depend on the choice of integration constants in the UV. Therefore its nu- 
merical value does not have an intrinsic meaning. However it can be used as a benchmark 
against which all other dimension-full quantities can be measured (provided one always 
uses the same UV boundary conditions). On the other hand, given a fixed set of initial 
conditions, asking that ?tio++ matches the physical value (in MeV) obtained on the lattice, 
fixes the value of ^ hence the energy unit. 

The holographic renormalization of such Einstein-dilaton theories is quite intricate as 
the AdS boundary conditions on the dilaton is unusual (0 — ?■ — oo near the boundary). It 
has been derived recently in []64| ]. 

3.1 Finite temperature 

In the large Nc limit, the canonical ensemble partition function of the model just described. 
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can be approximated by a sum over saddle points, each given by a classical solution of the 
Einstein-dilaton field equations: 



Z{f3) ~ e 



+ e 



-52 (/3) 



+ . . . 



(3.5) 



where Si are the euclidean actions evaluated on each classical solution with a fixed temper- 
ature T = 1/(3, i.e. with euclidean time compactified on a circle of length /3. There are two 
possible types of Euclidean solutions which preserve 3-dimensional rotational invariance. 
In conformal coordinates these are: 



1. Thermal gas solution, 



bl{r) [dr^ + df + dxmdx"") , $ = $o(r), 



(3.6) 



with r G (0, oo) for the values of P and Q we are using; 
2. Black-hole solutions, 



ds^ = b{ry 



fir) 



f{r)dr + dxradx" 



$ = <I)(r), 



(3.7) 



with r € (0, r/^), such that /(O) = 1, and /(r/^) = 0. 



In both cases Euclidean time is periodic with period (3o and /3 respectively for the thermal 
gas and black-hole solution, and 3-space is taken to be a torus with volume V30 and V3 
respectively, so that the black-hole mass and entropy are finite'^. 

The black holes are dual to a deconfined phase, since the string tension vanishes at the 
horizon, and the Polyakov loop has non-vanishing expectation value. On the other hand, 
the thermal gas background is confining. 

The thermodynamics of the deconfined phase is dual to the 5D black-hole thermody- 
namics. The free energy, defined as 



T = E-TS, 



(3.8) 



is identified with the black-hole on-shell action; as usual, the energy E and entropy S are 
identified with the black-hole mass, and one fourth of the horizon area in Planck units, 
respectively. 

The thermal gas and black-hole solutions with the same temperature differ at O(r^): 



6(r) = 6o(r) 



fir) 



1 



Cr^ 



where Q and C are constants with units of energy. As shown in [ 
enthalpy TS and the gluon condensate (trF^) : 

22 (tr F'^)t - (tr F^)^ 



+ ... r^O, (3.9) 
I they are related to the 



C 



TS 



3(47r)2 



240M3iV2 



(3.10) 



*The periods and 3-space volumes of the thermal gas solution are related to the black-hole solution 
values by requiring that the geometry of the two solutions are the same on the (regulated) boundary. See 
iH for details. 
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Although they appear as coefficients in the UV expansion, C and G are determined by 
regularity at the black-hole horizon. For T and S the relation is the usual one, 



T = -^, 5 = ^ = 4vr(MX^3)6^(r.). (3.11) 
For Q the relation with the horizon quantities is more complicated and cannot be put in a 



simple analytic form. However, as discussed in |38|, for each temperature there exist only 
specific values of Q (each corresponding to a different black hole) such that the horizon is 
regular. 

At any given temperature there can be one or more solutions: the thermal gas is always 
present, and there can be different black holes with the same temperature. The solution 
that dominates the partition function at a certain T is the one with smallest free energy. 
The free energy difference between the black hole and thermal gas was calculated in |38| 
to be: 



For a dilaton potential corresponding to a confining theory, like the one we will assume, 



the phase structure is the following |38| 



1. There exists a minimum temperature T^m below which the only solution is the 
thermal gas. 

2. Two branches of black holes ("large" and "small") appear for T > Tmm, but the 
ensemble is still dominated by the confined phase up to a temperature Tc > Tmin 

3. At T = Tc there is a first order phase transition to the large black-hole phase. The 
system remains in the black-hole (deconfined) phase for all T > Tc- 

The holographic mode has also been confronted successfully with recent lattice data 



|4C] at finite temperature, |39f|. 



4. Adding Flavor 

A number Nf of quark flavors can be included in our setup by adding space-time filling 
"flavor-branes" . In this case they are pairs of space-filling Z)4 — branes. 

To motivate the setup it is important to revisit the low-dimension operators (dimen- 
sion=3) in the flavor sector and their realization in string theory. At the spin-zero level we 
have the (complex) mass operator 

V^kVi ^ T,,- (4.1) 

dual to a complex scalar transforming as {Nf,Nf) under the flavor symmetry U{Nf)ji x 
U{Nf)L. At the spin-one level we have the two classically conserved currents 

They transform in the adjoint of the U{Nf)ji respectively the U{Nf)L symmetry. The 
flavor symmetry is expected to arise in string theory from Nf flavor branes (R) and Nj 



-16- 



flavor antibranes (L). Due to the quantum numbers, the vectors are the lowest modes of 
the fluctuations of the open strings with both ends on the D branes (^^), or the anti-D 
branes, Aj^. 

The bifundamental scalar T, on the other hand, is the lowest mode of the D — D strings, 
compatible with its quantum numbers. Its holographic dynamics is dual to the dynamics 
of the chiral condensate. This is precisely the scalar that in a brane-antibrane system in 
flat space is the tachyon whose dynamics has been studied profusely in string theory, pO| . 
It has been proposed that the non-linear DBI-like actions proposed by Sen and others are 
the proper setup in order to study the holographic dynamics of chiral symmetry breaking, 



1 48]. This dynamics was analyzed in a toy example, [55|, improving several aspects of the 
hard ||2^, and soft wall models, |^^. We will keep referring to T as the "tachyon", as it 
indeed corresponds to a relevant operator in the UV. 

The tachyon dynamics is captured holographically by the open string DBI+WZ action, 
which schematically reads, in the string frame, 

S[T,A'',A'^] = Sdbi + Swz (4.3) 
where the DBI action for the D — D pair is 



Sdbi = I drd'^x Str 



V{T) J- det {g^, + D{^TW,}T + FL^)+ (4.4) 



+ 



det {g^, + D^^TW^yT + pR)^ 



L R 

Here T is the tachyon, a complex Nf x Nf matrix. A^' are the world- volume gauge fields 
of the U{Nf)L X U{Nf)R fiavor symmetry, under which the tachyon is transforming as the 
{Nf,Nf ), a fact reflected in the presence of the covariant derivatives^ 

D^T = d^T-iTAj: + iA^T , D^T^ = d^T^ - iAj:T^ + iT^ A^ (4.5) 

transforming covariantly under 

T^VrTvI , A"- ^ VLiA^- - iVldVL)vl , A"^ ^ VRiA"^ - iVldVR)V^ (4.6) 

as well as the field strengths F^'^ = dA^^R — iAl^r A Al^r of the A^'^ gauge fields. 
A = e* = Nce'l' is as usual the 't Hooft coupling. We have also used the symmetric trace 
(= Str) prescription although higher order terms of the non-abelian DBI action are not 
known. It turns out that such a prescription is not relevant for the vacuum structure in 
the flavor sector (as determined by the classical solution of the tachyon) neither for the 
mass spectrum. The reason is that we may treat the light quark masses as equal to the 
first approximation and then in the vacuum, T = rl with r real, and this is insensitive 
to non-abelian ramifications. Expanding around this solution, the non-abelian ambiguities 
in the higher order terms do not enter at quadratic order. Therefore, for the spectrum we 
might as well replace Str — )• Tr. 



^We are using the conventions of | 
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The WZ action on the other hand is given by^'': 

Swz = Ti C A str exp [ilna'J'] (4.7) 

where M5 is the world-volume of the D4 -D4 branes that coincides with the full space-time. 
Here, C is a formal sum of the RR potentials C = i)^"^^^, and T is the curvature of 

a superconnection A. Note also that str in ([4.71) stands for supertrace and not symmetric 



trace. It acts on the space of D and D branes and is defined in Appendix C of [48|. 

In terms of the tachyon field matrix T and the gauge fields and A^ living respec- 
tively on the branes and antibranes, they are (We will set 27ra' = 1 and use the notation 

of m): 

\ T iAnJ ' \ DT iFr -TT^ ) ^ ^ 

The curvature of the superconnection is defined as: 

J='=dA-iAAA , dJ='-iAAT + iT AA = (4.9) 

Note that under (flavor) gauge transformation it transforms homogeneously 



In [^8[ the relevant definitions and properties of this supermatrix formalism can be found. 
By expanding we obtain 

Swz =Ti j C5AZo + C3AZ2 + CiAZ^ + C_i A Zq (4.11) 

where Z2n are appropriate forms coming from the expansion of the exponential of the 
superconnection. In particular, Zq = 0, signaling the global cancelation of 4-brane charge, 
which is equivalent to the cancelation of the gauge anomaly in QCD. Further, as was shown 



in ^ 



Z2 = dni , ni=istr{V{T^T))Tr{AL- AR)-logdet{T)d{StrV{T^T)) (4.12) 

This term provides the Stuckelberg mixing between Tr[Aj; - Aj}] and the QCD axion that 
is dual to C3. Unlike the 't Hooft limit, in the Veneziano limit this mixing happens at 
leading order in 1/Nc, Q. Dualizing the full action we obtain 



2iV2 



Scp-odd = ^J d^VaZiX) {da + i^if (4.13) 
j d'-x^ZiX) [d^a + xQd,y{T) - xV{T)A^f 



10 



This expression was proposed in and proved in BSl using boundary string field theory 
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with 

C = ^Imlogdetr , AL-An^^^AH + {Al-A%)\'' (4.14) 

and where we have set the tachyon to its vev T = rl . This term is invariant under the 
U{1)a transformations 

C^C + e , A^^A^-d^e , a^a-xeV{T) (4.15) 

reflecting the QCD U{1)a anomaly. It is this Stuckelberg term together with the kinetic 
term of the tachyon field that is responsible for the mixing between the QCD axion and 
the rj' . In terms of degrees of freedom, we have two scalars a, C and an (axial) vector, 
A^. We can use gauge invariance to remove the longitudinal components of A"^. Then an 
appropriate linear combination of the two scalars will become the glueball field while 
the other will be the rj' . The transverse (5D) vector will provide the tower of U{1)a_ vector 
mesons. 

The next term in the WZ expansion couples the baryon density to a one-form RR field 
Ci. There is no known operator expected to be dual to this bulk form. However its presence 
and coupling to baryon density can be understood as follows. Before decoupling the Nc 
branes, its dual form C2 couples to the U{1)b on the branes via the standard C2 A Fb 
WZ coupling. This is dual to a free field, the doubleton, living only at the boundary of the 
bulk. Once we add the probe L'4 + 1)4 branes the free field is now a linear combination 
of A^ and an Nf/Nc admixture of A^ originating on the fiavor branes. The orthogonal 
combination is the baryon number current on the fiavor branes and it naturally couples to 
Ci. Therefore the Ci field is expected to be dual to the topological baryon current at the 
boundary. 

Finally the form of the last term requires some explanation. By writing Zq = dQ^ we 
may rewrite this term as 

'foAQ^ , Fo = dC_i (4.16) 



Fq ~ Nc is nothing else but the dual of the five-form field strength. This term then provides 
the correct Chern-Simons form that reproduces the fiavor anomalies of QCD. Its explicit 
form in terms of the gauge fields ^l,_r and the tachyon was given in equation (3.13) in p8| . 

The action as described is based on the fiat space Sen action for the D — D brane- 
antibrane system. In the presence of curvature and other non-trivial background fields, 
like the dilaton we expect corrections to the DBI action. Such corrections may affect the 
tachyon potential as well as the kinetic terms of the vectors and the tachyon. Some generic 
properties are expected to remain though, and these include the tachyonic nature of the 
scalar near the AdS boundary and the exponential asymptotics of the potential at large T. 

5. The bottom-up models 



At X — -jj- — the IHQCD model, [|36|] , is described by the action 



Sg = M^Nl j d^x^ 



R-l^^ + V,iX) (5.1) 
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with A the 't Hooft couphng (exponential of the dilaton) and a potential that has the 
following asymptotics. 

limVgiX) = — [l+viX + V2X'^ + ■■■] , lim yg(A) ~ Atyi^ (5.2) 

At finite x we must add the fiavor sector. For the vacuum structure, and with all 
masses of the quarks being equal it is enough to add the C/(l) part of the tachyon DBI 
action, 

Sf = -xM^N^ j d^x Vf{\,T)yJdei{g^^ + h{X,T)d^TduT^) (5.3) 

where we have set the gauge fields to zero. The total action is 5" = 5"^ + Sf. Note that the 
overall sign of the DBI action is negative. The function Vg and its asymptotics has been 
discussed in detail in |3f:]. We will consider it known, and when needed we will use the 
form that was in agreement with YM data, |39]. The tachyon potential Vf{X,T) should 
satisfy some basic principles. For flat space D-branes, Vg ~ je~^^'^^ . In our case, near the 
boundary, where T — 0, A — >• 0, we expect, in analogy with Vg, a regular series expansion 
in A, T 

Vf ~ VoiX) + Vi{X)T^ + O(r^) (5.4) 

with Vo,i(A) having regular power series expansions in A. As we will see later, the functions 
V(),i(A) may be mapped into the perturbative /3-functions for the gauge coupling constant, 
and the anomalous dimension of the quark mass operator. 

Near the condensation point, T — )• cx) we expect the potential Vf to vanish exponen- 
tially. This is based on very general arguments due to Sen that guarantee that the brane 
gauge fields disappear beyond that point. 

Finally the function /i(A, T) was introduced to accommodate the fact that the action 
in (5.3) is written in the Einstein frame. In fiat space, this factor is unity in the string 

4 

frame but becomes nontrivial (/i ~ A~3) in the Einstein frame. 
5.1 The equations of motion 

Collecting the action of the glue and fiavor sectors together, 

^2 



C = (M^ 



^{r- + VgiX)^ - X V>(A, T) Vdet {gab + h{X,T)daT d^T) . 

(5.5) 

We shall take the following Lorentz-invariant Ansatz for the metric: 



ds^ = e^^ {dxl^ + dr^) . (5.6) 

In our Ansatz the warp factor A, the scalar A and the tachyon T depend only on the radial 
coordinate r. 

The Einstein equations take the form: 



Rab-\9abR = Tl, + Ti,, (5.7) 
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where T^^ and T^^ are the energy momentum tensors of the glue and flavor sectors, respec- 
tively. These equations translate into 



12(A 



i\2 



,2A 



e-2^/i(A,r) {T'f 



where primes stand for r-derivatives. For x = these equations agree with 
the equations of motion for the dilaton and the tachyon are given by: 



A" - + 2, A X' 



■ a2 e^^ 
X dh 



dVn 



f 



dX dX 



l + e-2^/i(A,r) {T 



'/\2 



(5.8) 

(5.9) 
Finally, 

(5.10) 



2 dX ^l + e-2^/i(A,r) {T'y 



r^ii OA f .1 hX! X' dh\ 

/ A' 5Ff X' dh\ . 

6. Conformal fixed-point solutions 



1 dh 1 dVi 



2hdT Vf dT 
e^^ dVf 



{T'Y+ (5.11) 



h Vf dT 



0. 



There are solutions to the equations above that are conformal, and are related to the fixed 
points of an "effective potential" . To find them we must set X' = T' = 0, and A" = T" = 
in the equations (5.10), ( ^.11 ) which imply that we must be at a critical point of the 
effective potential. 



V,s = Vg{X)-xVf{X,T) , drV^s = dxV,s = 
For any solution A*,T=k of the above conditions, we obtain an AdSs space with 

72 = K!ff(A*,r*) 



(6.1) 



(6.2) 



as is obvious from equations (5^) and ( |5.9D . 

When Vf has the standard dependence on the tachyon [Q, there are two solutions to 
the condition dxVes = 0, namely T = (chiral symmetry unbroken) and T — )■ oo, (chiral 
symmetry broken). 



T = 0. In this case the second condition of extremality ( |6.1[ ) is dx{Vg{X)—xVf{X, 0)) = 
0. In the UV, A — )• 0, this potential is constructed to emulate the perturbative (3- 
function, and therefore has a free-field theory fixed point. This is the UV theory. In 
the IR, it will also have a fixed point at A = A,,,. In the BZ region A,,, ^ 1. A priori 
there are two possibilities. 



1. The fixed point disappears for x < x^, for a given x^,. 



-21- 



2. The fixed point exists for all < x < -j-. 

We will discuss these two options later on in this paper. 

• T — )• oo. In this case as Huit-^oo ^(^i T) = we obtain that the second extremality 
condition is dxVg = 0. This is equivalent to the existence of a fixed point in large-N 
YM theory. It has been argued however in that this is not true. Therefore, there 
is no fixed point with T — t- cxd. 



7. Holographic /3-functions 

In holographic theories we may define non-perturbative /3-functions that capture the de- 
pendence of coupling constants on the RG scale. This concept was developed and used first 
in [^], in order to explore the physics of Einstein dilaton theories and their holographic 
relation to YM theory. In particular in pH] it was shown that the /3 function is intimately 



related to the generalized superpotential. It was shown later in |38| that such defined 



/3-functions are indeed related to the quantum breaking of scale invariance, and appear in 



the trace of the stress-tensor. Such relations were shown in full generality in |64| and have 



been confirmed also in |65|. 



To define the /3-functions we need a notion of energy scale. At the two-derivative 
level, such a function is the scale factor e^, and indeed near the AdS boundary it can be 
identified as the energy scale. It remains always a decreasing function, and becomes zero 
in the ultimate IR. It can therefore be taken as the energy scale in the whole of the bulk 
space. 

We therefore define the /3- function and "anomalous" dimension 7 as^^ 

^.,(XT) , ^.,(A,T) (7.1) 

The equations of motion provide equations for the /3-functions. To obtain them we must 
convert radial derivatives to derivatives with respect to A, 

X' = A'/3 , T' = A'j , X" = i3A" + {j/3T + f3f3x)A'^ , T" = + (77T + /37a)^" 

(7.2) 

where /3t = ^ etc. Substituting in (^^) we obtain 



^^Notice that, as we shall see later, it is the ratio ^{\,T)/T (rather than 7(A,T)) which corresponds 
closely to the anomalous dimension of the quark mass in QCD. Excluding the extra T in the definition 
simplifies many of the equations below. 
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which can in principle be solved for as a function of A, V^, Vj, /i, /3, 7. In general, the 
solution is not unique. Similarly from ( |5.8[) we obtain 



72A 



1 + 



9A2 



p2A 



+ 



2/32 W'2 Vg-xVfJl + h^^^ 



Then ( 5.10 ) and ( 5.11 ) become 



(7.4) 



a/3 Bp /32\ A 



1/2 



^^5A A j ^ 6 



(7.5) 



^A2 



1 91og_t^ / 



aK ay, 



, a; a/i yf72e-2^^'2 



^ 2" aA 



1/4 



^7' -IT (7-6) 



+ 



/3'V.^log(^\ /I 



(9 log /i (9 log V/ \ 2 



ar 



ar 



7 



I2A 



where we eliminated j4" by using ( [7. 41 ) . This is a system of two first-order partial differential 
equations for /3, 7, with inputs Vf,Vg,h. The equations are highly non-linear. Setting x = 0, 
the system reduces to 

/32 \ 



12 1 



9A2 



^ar ^^aA a / 8 dx 3 ^ 



In this case ^ = and ( |7.7D can be rewritten as 



9A ~ A 



4 + 



DA^aiog^ 
2/3 dX 



1 



9A2 j 



(7.7) 
(7.8) 

(7.9) 



These equations are similar to those in the probe (quenched) limit. 

Note also that although the equations of motion are linear in x, the /3-system ( [7^5| ,[fy 
is non-linear in x, with the understanding that 4i is to be eliminated using (|7.3|). 



7.1 The UV Fixed point 

Near the boundary where A, T — t- 0, we expect that 7 — ?■ 0, so that equation ( |7.3| ) becomes 



/32 \ ^'2 



A'^ _ Vg-xVfiX,0) 



p2A 



12 1 



/32 

9P 



+ 



(7.10) 



^Standard "linear" scalar theories with action given by a kinetic term plus a potential have a unique 
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According to the discussion of Sec. |3|, the potentials are expected to be regular at A = 0, 

Vg = Vq+ViX+V2X^+0{X^) , xVf = Wo+WiX+W2X^+Oi\^) + {Zo+ZiX+Z2X'^)T^+- ■ ■ 

(7.11) 

h = ho + hiX + h2X^ + O(A^) + 0{T^) 
We also take the following Ansatz for the beta functions: 

f3 = -boX'^ + biX^ + --- + OiT^) + --- , 7 = (7o + 7iA + 72A2)r+--- (7.12) 

Inserting these into the equations (|7.5| , [7.6D , we find 



_9 Vi-Wi ^ 207 {V2 - W2f 9 Vi-Wi 

° SFo-Wo ' ^ 256 (Fo-VFo)2 ^.Vq-Wq ^ ' ' 



12Zo ( Zi Wi hi Vi-Wi 

71 



(70 + 2)hMVo - Wo) V Zo ho Vo - Wo 

To leading order the solutions around the UV fixed point are 

\ = ^ + boA + --- , r = Toe^°^ (7.15) 

In order for T to have the proper UV dimension we must have 70 = —1. In the massless 
case, T is dominated by the vev, and we have to choose 70 = —3. These solve Eq. ( 7.14 ) if 

2iZo 



hoWoiVo - Wo) 



-3 . (7.16) 



As we shall point out below, the combination 7(A, T)/T is mapped to the anomalous 
dimension of the quark mass in QCD ( |2.8D . Remarkably, the solution ( 7.12| ) is consistent 
with QCD perturbation theory: 'y{X,T)/T has a series expansion in A. The leading term 
is fixed according to the UV dimension, and the correction terms are identified with the 
anomalous dimension. 

7.2 Confining IR asymptotics 

In the IR, a confining asymptotic has the property that ^ — )• 00. There are two possibil- 
ities for the tachyon: 

1. If T = 0, 7 = and then Vg Vg - xVf{T = 0). In this case 

^2 ^ Vg-xVf{X,Q) 



+ ••• (7.17) 



2. T — )• cxD in the IR, — )• exponentially, and we obtain the pure YM case in the 
IR. Here 
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If we assume that Vj(A,0) does not grow faster than Vg{X) in the IR, the solution in both 
cases is the same as in pure YM (x = 0), 



f3 



3A 



3 1 

^ 8 log A ^ 



(7.19) 



This non-perturbative /3-function indicates that the scahng dimension of the dual operator 
tr[F'^) in the IR is A = | read from its linear term. However this does not imply there is 
a scaling regime in the IR, as the metric is far from AdS. 
The tachyon equation becomes 



d log Vj 
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dT 



+ h[ Vg-XVf\ 1+h-f 
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+ 4 + /3 



dlogVf pdlogh\ h^V^ 



dX 
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2 dX 
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dT 



+ /3^+(2(l 



9A2 



+ /3 



diog{hVf) \ 



dX 



1 dlogh d logVf , 2 
7+ i;^^ 7^;^ 7 



/32 

(7.20) 



2 dT 



dT 



hVa 



12 1 



/32 
9A^ 



In the IR, for a class of potentials 7 — )• —00 and a/i + /i7'^4^ 



7. We ex- 



pect that Vf\ 1 + /i7^^^ — )• as V/ vanishes exponentially while 7 increases polynomially 



with increasing T. We also expect hVg to be approximately constant and ( 1 — ^ 



We also expect ^^^^ — )• and to be subleading. 

Therefore the leading terms in equation (7.20) are expected to be 

d \ogVf p diogh X hVg 73 dj_ _ d logF/ 2 ^ 
5A 2 dx J 9 ^ar ^ 

If we define 6 = 6 satisfies a linear equation 
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4- 



(7.21) 



d6_ dlogVf _ hVg_ 
df dT ~ ~9~ 



4 + /3 



d iog{VfVh) \ 

dX J 



(7.22) 



with solution 



C hVg 



with T* large enough so that we are in the IR regime. For 7 we obtain 



7 



C+!^JlVf{A + l3^^2^^)dT 



(7.23) 



(7.24) 



Note that if /i ~ A"!, then /3 ^isg^ ~ l. If 

we also assume that Vf has factorized 

dependence, Vf = Vfx{X)VfT{T) we obtain 
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Note that from the whole one-family of solutions above, only one is diverging for large T 
and this is the only reliable solution. 

For a tachyon potential of the form VfT = e~"''^^ the diverging solution is 

T + .... (7.26) 



For the type of potentials we will be using hVg / yjlog A approaches a constant which we 
denote by b, and Vf\ ~ so that 

18a ^ 9a ^ , 
7Ci , ^r + ---~ l=T^ (7.27) 

Note that if we assume a more elaborate tachyon potential of the form Vj = yyo(A)e~"^'^^^^ 
with a (A) increasing with A in the IR then the asymptotic behavior of 7 changes and it 
vanishes in the IR. To capture this behavior from ( 7.20| ) we take 7 — )• in the IR and 



1 + /i7^^^ — 7" 1. We also again expect that hVg is approximately constant, ^1 — ^ 

|, and ^g^^ — )• 0. We also expect that the derivative terms 7^ and are subleading. 
The two leading terms in ( [7. 20] ) are then expected to be 

9 logy; hVg p d\og{hVf) _ 

— WF- + —9 d)^^ - ° ^^-^^^ 

The solution behaves for large T as 

18a(A) 1 12 1 

'^"/3W,a'(A)r + ---"-6|M^^r + --- ^ 

where h is defined as above. 
7.3 Some simple /? functions 

We may engineer a (3 function that interpolates between the one-loop perturbative YM 



/3-function, and the non-perturbative one in ( 7.19 ). We could also have a 7 function 



interpolating between an operator with dimension Auv in the UV, and Air in the IR: 

R(x\ u ,,^ l + (AiR-4)r2 

/5(A) = -^0^^ ' ^ = (^^^-'^^ l + (Aw-4)r^ ^^-'"^ 

The flow equations can be integrated to 

---logA = 6o^ , ^^-—^^ + ^^^^^ _ ^^^^^^ _ log (1 + (Air - ,)T ) = A 

(7.31) 

Such /3 functions can be converted into potentials Vg,Vf, via the equations ( |7.5D , ( |7.6|) . 



-26- 



Figure 1: The solutions for tlie (3 (top row) and 7 (bottom row) functions ioY x — 2 (left), x — 3.9 
(middle), and x = A.2 (right). We added the red lines on the top row at /3 = in order to show the 
location of the fixed point. See the text for a detailed explanation. 



7.4 Numerical solutions for the (3 and 7 functions 



The equations ( |7.3| ), ( [7.5D , and ( |7.6D for the /3 and 7 functions can be solved numerically 
for fixed potentials by combining solutions evaluated along the RG flow, as detailed in 
Appendix Fig. ^ shows the results for various values of x. We used the potentials of 
scenario I from Appendix |C| and required that the solutions flow to the good IR singularity, 
as explained in Appendix ^ Notice that the values x = 3.9 and 4.2 were chosen to lie 
slightly above and below the edge of the conformal window, which is at Xc — 3.9959 for 
the potentials used in the plots. 

The plots of the beta functions show a smooth transition as T evolves from small 
^ 1 to large > 1 values, reflecting the expectation of Sec. |6|. For small tachyon, the beta 
function has a fixed point corresponding to the maximum of the effective potential Ves, 
which moves to lower values of A as x is increased, whereas for large T the A-dependence 
of the beta function approaches the Yang-Mills form. The 7 functions show a transition 
between the small and large T regions as well. For large T and A, the solutions agree 
with Eq. ( [7.27 ). For small T, the structure is richer. For x = 2 the 7 tends to a constant 



value as T — )• (except for very small A) so that "f{X,T)/T, which is plotted in Fig. |T|, 
diverges. This behavior is pushed for larger A as we increase x to 3.9, and has disappeared 
for X = 4.2, and the gamma function is instead linear in T in accordance with Eq. ( 7.12| ). 
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These changes are related to the tachyon changing sign during the flow. If the tachyon 
has a zero, we hit the T = line before reaching the UV singularity. Then 7 approaches 
a constant value as T — 0, and Eq. ( 7.12 ) does not apply. As will be discussed in detail 
below, tachyon zeroes are indeed expected for x < Xc- Notice also that j/T tends to —1 as 
A —7- with fixed T at least for small T, reflecting the expected value of 70 in Eq. (7.12). 



8. The background solutions 

8.1 Generic properties of the background 

We start by discussing the symmetries and integration constants of the equations of motion 
of V-QCD, (p|)-(l5Tl|). We shah assume the exponential Ansatz 

Vf{X, T) = Vf^{\) exp [-a(A)r2] (8.1) 

for the potential of the tachyon DBI action. Then the background EoMs have the following 
symmetries 

1. A ^ ^ + log A , Vg^ K-^Vg , Vfo ^ A- V/0 (8.2) 
T ^ AT , a ^ A^^a ; 

2. T-^AT, a^A^^a, h-^A^'^h; (8.3) 

3. r A(r - ro) , A-logA . (8.4) 

The first symmetry can be used to fix the value of the UV AdS radius i, and is usually 
associated with the units of energy in the boundary theory. The second one will be used 
to fix the normalization of h in the UV. The third one is essentially different from the first 
two, since it does not involve the potentials. It will therefore remain as a true symmetry 
of the background solutions. 

We may choose a set of independent equations of motion which contains one first order 
and two second order differential equations. Therefore, their general solution includes five 
integration constants. These can be identified as the coefficients of the UV expansions 
of the fields as follows (assuming that the solution has the standard UV singularity with 
A — >• 0, see Appendix |D[) . The tachyon UV expansion has the usual free constants related 
to the normalizable and non-normalizable solutions, identified as the quark mass m and 
the vacuum expectation value a of the qq operator, respectively. 

In close analogy, the solution for A involves two constants, identified as the UV scale 
A = Auv of the expansions, and another constant A which we will define in Section |9|, 
related to the gluon condensate and the free energy of the system. The fifth integration 
constant is simply the location of the UV singularity which can take to be r = by the 
translation symmetry of Eq. ( ^.4|) . We shall require that the system has a repulsive, "good" 
kind of IR singularity, which fixes the values of the condensates a and A in terms of m and 
Auv- 

In addition, we still have the scaling symmetry of Eq. (|8.4| ), which can be used to vary 
the units of all constants, and refiects the corresponding scale transformation of the dual 



-28- 



field theory. Therefore, the single parameter which characterizes all nontrivially linked 
physical backgrounds (for fixed potentials and, in particular fixed x = Nf/Nc) is the ratio 
of the "source" coefficients m/Auv^^- 

For the choices of the functions h, Vg, and Vj of interest to us, the tachyon typically 
decouples from the other fields asymptotically in the UV and in the IR. The UV and IR 
asymptotics are discussed in detail in Appendices ^ and 0, and we shall repeat only the 
main features of the physically interesting possibilities here. The physically relevant UV 
asymptotics are restricted. As pointed out in |36|, the fields A and A can be expanded in 



— 1/logr at the UV boundary r = in the probe limit. 

Similar expansions work also at finite x = Nf/Nc- The tachyon is required to vanish 
linearly T(r) ~ mr (or faster if m = 0) in r in the UV. Taking e = —1/logr — )• 0, the 
tachyon T ~ mexp(— 1/e) vanishes exponentially while A and A have power-like behavior 



in e. Since, in addition, the functions h and Vf must be regular in the UV (see Sec. |8.2 
below), the tachyon can be set to zero in the leading order action. We find that A and A 
satisfy their probe limit equations of motion, but with the dilaton potential Vg replaced 
by Uefr(A) = V^(A) — xV/(A, T = 0), which also verifies that the UV expansions of A and A 
have the same form as in the probe limit. 

We shall only discuss cases where the tachyon indeed decouples asymptotically in the 
IR. This is the case, if we take a tachyon potential having an exponential T dependence, 
V/(A,T) DC exp(— aT^), and the tachyon has a power-law (or faster) divergence as r — )■ oo. 
This is indeed what is suggested by tachyon condensation in string theory. Therefore, the 
flavor part of the action is exponentially suppressed in the IR. Consequently, the tachyon 
decouples asymptotically from A and A, and their asymptotic expansions in the IR have 
exactly the same form as in the probe limit, and are determined by the potential V^(A). 

In summary, even though all fields couple nontrivially for general values of the co- 
ordinate, the probe limit description will be valid in the UV and IR. In particular, this 
guarantees that the interpretation of the integration constants is the same for finite x as 
in the probe limit. The decoupling of the tachyon in the IR leads to the system having 
similar "good" IR singularities as in the probe case. 

An important difference with respect to the probe limit discussion is that the potentials 
which characterize the backgrounds in the UV and IR regions will be, in general, qualita- 
tively different. As discussed above in Section the potential V^fr(A) = Vg{X) — xVf{X, T = 
0) which controls the UV behavior will be chosen such that it admits a fixed point (ex- 
tremum of the potential) at least for large values of x as required by the Banks-Zaks 
analysis. 

The fixed point of V^s will play an important role in the dynamics in the intermediate 
region between UV and IR. In particular, for identically vanishing tachyon, the background 
simply fiows from the UV fixed point at A = to the IR fixed point at finite A. Adding, a 
tiny quark mass (or chiral condensate), the solution in the UV region will not be changed 
drastically. However, no matter how small the quark mass is, the tachyon will eventually 



^^Due to practical reasons we shall often use a parameter To linked to the tachyon behavior in the IR 
instead of the quark mass to characterize different backgrounds. 
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become large and start coupling to A and A, which will drive the system away from the IR 
fixed point. 

We can make two observations. First, the special case of zero tachyon will be essen- 
tially different from all other solutions. This solution has vanishing quark mass and chiral 
condensate, and is therefore identified as the solution corresponding to conserved chiral 
symmetry. We will discuss this special case separately below. 

Second, there is a possibility to have solutions which come very close to the IR fixed 
point, but are eventually driven away by the increasing tachyon. Such solutions will be 
identified with quasiconformal or "walking" dynamics of the dual field theory, where the 
coupling constant remains approximately fixed over a large range of energies. We shall see 
below how the phase structure of QCD in the Veneziano limit, which includes a quasicon- 
formal region, arises in our class of models. 

8.2 Matching UV behavior w^ith the QCD /3-functions 

We shall now discuss the most important links of the potential functions to the physics of 
the dual field theory. We start by an analysis of the UV region, where the behavior of the 
system can be mapped to the QCD /3-functions [^] as already discussed in Section |^. In 
particular, in the probe limit (a; — )• limit), the UV behavior is controlled by Vg{X), that 
has the expansion 

Vg{X) = Vo + ViX + V2X'^ + --- (8.5) 

Here Vq > can be freely chosen and it fixes the AdS scale for x = as Vq = 12/£q. The 
other coefficients can be mapped to the Yang-Mills /3-function. At one-loop order we have 

l3e 



12 

^0 



i.6) 



where bo^^ is the one-loop coefficient of the /3-function, from which Vi can be solved. 

Moreover, as discussed above, at finite x the UV behavior is similar to the probe limit, 
but the role of Vg is taken by the potential V^g (A) = Vg{X) — xVf{X, T = 0). Therefore, we 
take 

Vf{X,T = 0) = Wo + WiX + W2X^ + --- (8.7) 
and the relation to the QCD /3-function reads 



12 

V,siX) = VgiX)-xVf{X,0) = - 



1 + -boX + 
y 



where bo is the leading coefficient of the /3-function in the Veneziano limit. 

Similarly to Vq, the coefficient Wq can be freely chosen, and the other coefficients can 



be solved from Eq. B.8. However, there are constraints: the AdS scale must remain positive 
for all < X < 11/2, and Wq should also be positive (see Appendix j^). These boil down 
to 

< H^o < ^Vo . (8.9) 
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In addition, as pointed out in Section |7[ we can map the UV behavior of the tachyon 
action to the (anomalous) dimensions of the quark mass and the chiral condensate of the 
dual field theory. For definiteness, we parametrize V/(A,T) = V/o(A) exp(— a(A)T^) and 
assume that h depends only on A. Then the dimension of the quark mass constrains the 
UV behavior of h{X) and a(A). A detailed analysis is done in Appendix p.l.l| . 

Remarkably, assuming that the potential functions have analytic expansions at A = 
is consistent with QCD perturbation theory. Indeed, requiring the dimension of the quark 
mass to approach one in the deep UV fixes the leading terms as 

h{\) 



a{X) 



2P 

— {l + hiX + 



10) 



and the next-to-leading coefficient hi can be matched with the one-loop anomalous dimen- 
sion of the quark mass 7™, (A). 

By using Eq. ( D.1C| ) from Appendix, we obtain 



9 

= o 



48, 4^ ■ 
3 9^° + 3^^ 



Ul) 



where 70 is the leading coefficient of the anomalous dimension, — dlog m/dlog = 7m(A) 
7oA+- • • . Notice that, for example, the non-normalizable term in the tachyon solution (| 
reads after this identification 



T{r) = mr(-log(rA))-^»/''« 



1 + 



1 



log(rA) 



^.12) 



where the logarithmic correction is consistent with the one-loop solution for the running 
quark mass in QCD. 

8.3 Condensate dimension and the edge of the conformal window 

The most important new feature of the system discussed in this article is the description 
of the phase diagram of QCD as a function of x = Nf/Nc in the Veneziano limit. We 
shall now indicate which potentials lead to the desired structure. This constraint to a large 
extent independent of the one discussed above that was set by the UV expansions. 

The phase structure is linked to the dimension of the chiral condensate at the IR fixed 
point which is found at the maximum of V^s = ~ ^^fo in the conformal window (in the 
limit of small tachyon background). 

We denote the value of the coupling at the fixed point by A* so that 



v;'fr(A) = y;(A,)-xy;o(A*) = o. 

From the action (^.5[) we can calculate the IR AdS scale 



12 

w 



(8.13) 



(8.14) 



To calculate the dimension we need the mass of T in the IR fixed point, which can be 

IR ~ ~~r7V~r ' \o.l5) 



extracted from the action (5.3) and reads 



m 
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where we again parametrized Vf{X,T) = Vjo(A) exp(— a(A)T^). From this we obtain 



Air (4 - Air) 



-miR^^R 



24a(A,) 



.16) 



The desired phase diagram is only obtained if the expression Air (4 — Air) 



-m 



has a certain dependence on x. It must start at a value smaller than 4 at the end of the BZ 
region (which is guaranteed as the standard UV boundary conditions fix it to 3 there), grow 
as X decreases, and become 4 at some value of x so that the BF bound is saturated 
0. Indeed, the solutions near the edge of the conformal window stabilize such that the 
critical Xc, defining the location of the conformal phase transition for massless quarks, is 
determined by 



Air (4 - A 



IR, 



^.17) 



We shall discuss why this is the case in detail later on. Note that in addition to the 
explicit dependence on x in the factor xVfo, Eq. (^.16 ) depends on x through A* due to 
the definition of Eq. ( 8.13| ). 



We stress that the saturation of the BF bound means that for theories near the critical 
Xc, Air — )• 2, or in other words, the anomalous dimension of the quark mass at the fixed 
point — )• 1. That is, our model reproduces the standard assumption for the energy 
dependence of the chiral condensate near the edge of the conformal window. Recall that 
this is extremely important for realizations of walking technicolor. 

If the potentials are matched with the UV physics of QCD, the expression ( |8.16| ) can be 
further simplified. Fixing the tachyon mass in the deep UV we obtain a(0)//i(0) = 3/(2£^). 
In the limit T — )• 0, the beta function /3(A) = dX/dA satisfies a first order differential 
equation which depends on Veff{X), as can be immediately concluded by comparing our 
action to the probe limit one [^6| for T = 0. In terms of the phase function X, defined by 



we find 



A 



^ = /3(A) = 3AX(A) , 

dxix)_ fsxw + sx'-^^^^)'-^^ 



dX \^--v-v — • J gj^(^^ 

Solving this, we can parametrize the potential in terms of the /3-function: 



V,s(.X) = Vg{X)-xVfo{X) = ^ 
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The IR dimensions now satisfy 
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Figure 2: The squared tachyon mass at the IR fixed point Air (4 — Air) = —mf^ii^ as a function 
of X for the two scenarios described in Appendix ^ Thick blue, and thin red curves are the results 
for the scenarios I and II of Appendix |c], respectively. The dashed lines show the maximal changes 
as Wo is varied from (upper curves) to 24/11 (lower curves). 



where in the second line we used the fact that /3(A*) = 0. If we do the matching with 
perturbative QCD results in the UV, the /3-function can be identified as the QCD /3- 
function in the Veneziano Hmit. Then the exponential factor is small in the sense that it is 
roughly proportional to boX* where both the one-loop coefficient of the /3-function bo and 
the coupling at the fixed point A* vanish in the BZ limit x — )■ 11/2. Therefore we expect 
that Air (4 — Air) depends on x mostly through the functions h and a. In the two- loop 
approximation of the /3-function we obtain 

^.„A. . (^22) 



A2 9 

In Fig. § we plot the dependence of Air(4 — Air) on x for the two scenarios of potential 
choices of Appendix]^. The solid lines give the results for the value Wq = 12/11 used in 
the Appendix, and the dashed lines show the sensitivity of the result for the choice of Wq 
as this parameter is varied over its allowed range. In all cases. Air (4 — Air) intersects the 
value of 4, shown as the horizontal dotted line in the plot. This suggests that the class of 
potentials that has the desired phase structure quite in general includes those ones that are 
matched with the UV behavior of QCD. Moreover, the critical value of x is found within 
a quite narrow band 

3.7 < X, <4.2 , (8.23) 

and the largest source of uncertainty is the choice of Wq-^^ We stress that there is no strict 
bound on the allowed values of Xc, but the numbers in Eq. ( |8.23 ) rather give the expected 



magnitude for the variation of Xc within natural potential choices. 

It is also instructive to show the dependence of the anomalous dimension of the quark 
mass at the IR fixed point on x (see Fig. |3|). Within our model the anomalous dimension 
is defined by 7=,, = Air — 1 where Air is the smaller of the two roots. The result for the 



*We have checked this also for some additional potentials that are not discussed in this article. 



-33- 



r* 



1.0 




0.8 



0.6 



0.4 



0.2 



X 



4.0 



4.5 



5.0 



5.5 



Figure 3: The anomalous dimension of the quark mass at the IR fixed point as a function of 
X within the conformal window in various approaches. The solid blue curve is our result for the 
scenario I of Appendix The dashed blue lines show the maximal change as Wq is varied from 
(upper curve) to 24/11 (lower curve). The dotted red curve is the result from a Dyson-Schwinger 
analysis, the dot-dashed magenta curve is the prediction of two-loop perturbative QCD, and the 
long-dashed green curve is based on an all-orders /3-function. 

potentials of scenario I of Appendix ^ is shown as the solid blue line in Fig. ^ and as in 
Fig. |2|, the dashed lines show the maximal variation as the parameter Wq is varied over its 
allowed range. 

Our result is very similar to the prediction obtained by calculating the anomalous 
dimension from Dyson-Schwinger equations in the rainbow approximation |^^, and evalu- 
ating the result at the zero of the perturbative two-loop /3-function (|2.3D (dotted red curve 
of Fig. ^). As in our model, the conformal window ends at the point where 7^, reaches one 
and becomes complex in this approach. The deviation from the simple perturbative esti- 
mate (dot-dashed magenta curve), which was obtained by using the two- loop anomalous 
dimension (|2.8D instead, is also small. However, an all-orders /3-function []6^ ] (long-dashed 
green curve) predicts much smaller values of 7* for low values of x. 



8.4 Constructing the background solutions 

We will select concrete potentials for evaluating the background numerically. These po- 
tentials must satisfy the constraints of the previous two subsections in order to produce 
the desired phase structure of QCD in the Veneziano limit as well as the perturbative 
UV physics. In addition, the system needs to have an acceptable IR singularity where 
the tachyon diverges, which adds extra requirements to the large A behavior of the poten- 
tials. However, as we shall demonstrate briefly, different choices for the IR behavior do not 
change the qualitative features of the background. 

The construction of explicit potentials is detailed in Appendix P, where also the IR 
behavior is fixed by using the generic analysis of Appendix ^. For clarity, we repeat the 
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final result here. The potentials of the scenario I in Appendix |C| read 
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We shall use this choice when calculating the background numerically, unless stated oth- 



erwise. 



Recall that only three of the equations of motion (|5.8|)-( |5yTl ) are independent. We 
choose a set of three second-order equations for the numerical calculations, and treat the 
remaining first order equation as an extra constraint. As usual, we choose boundary 
conditions that satisfy the constraint, which is then automatically satisfied for all r}^ 

When solving the equations numerically, it turns out that shooting from the IR is nu- 
merically stable. We shall first discuss backgrounds where the tachyon has a nontrivial 
profile. We fix the boundary conditions in the deep IR by using the asymptotic IR expan- 
sions at the "good" IR singularity at r = oo (see Appendix E.2.2). For the functions given 
above, the asymptotics becomes 
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Here we already used the translation symmetry of Eqs. ( ^.4| ) to set ro in the formulas ( E.l^ ) 
and (F^ ) to zero, and we can further use the scaling symmetry to fix R= 1. After this. 



^^In all numerical calculations, we shall do the coordinate transformation from r to A discussed in 
appendix ^ because after this, the UV structure of the solutions is reproduced more accurately. However, 
as this transformation is straightforward, we shall continue to discuss the solution in the "r-space" . 

^® There is a potential instability related to the constraint, as it is exactly satisfied only near the IR cutoff 
due to numerical effects. For general r there is an error which may grow exponentially as we solve the 
system towards the UV. However, we have the freedom of modifying the system of second order equations, 
which is used to calculate the solution, by adding multiples of the constraint to some of the equations. In 
this way the error can be made to decrease exponentially instead of growing towards the UV, so that the 
instability is removed. 

^^Notice that with this choice of fixing the translation invariance, the UV boundary of the final solution 
will not be at r = 0. After obtaining the solution, we can relax this condition and use translation invariance 
again to move the UV singularity at r = 0, if desired. 
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Figure 4: The UV behavior of the background solutions with good IR singularity for the scenario 
I (left) and scenario II (right) defined in Appendix |c[ The thick blue curve represents a change in 
the UV behavior, the red dashed curve has zero quark mass, and the contours give the quark mass. 
The black dot where the zero mass curve terminates lies at the critical value x = Xc- For scenario 
I (II) we have Xc — 3.9959 {xc — 4.0797). See the text for detailed explanation. 



the solution (for fixed x) depends nontrivially on only one free parameter Tq in the IR, as 
is characteristic for a good IR singularity. 

To continue, we choose an IR cutoff -Rir such that the tachyon is large, and therefore 
decoupled from the fields A and A. With the above choices of potentials, a sufficient value 
turns out to be r(i?iR) = 70: with this choice, the tachyon has decoupled to a good 
precision, and Riyi/R is large enough for the asymptotic expansions to work. We use the 
asymptotic expansions to fix the values of A, A, T, A' , and T' at r = Rib., and solve A' 
from the constraint (the first order EoM). The solution is then obtained by numerically 
solving the set of second order equations of motion toward the UV, until some of the UV 
asymptotics described in Appendix ^ are reached. 

The obtained UV behavior is depicted in Fig. ^ (left) as a function of the only remaining 
free parameters, Tq and x. For comparison, we also present the same plot for the scenario 
II of Appendix on the right hand side, where Tq is replaced by the parameter ri of 
Eq. ( p.21 ). For clarity, we shall only refer to the variable Tq in the discussion below. 
Because of the invariance of the Lagrangian under T — )• —T, scanning over positive Tq is 
enough to catalog all possible solutions. We shall explain the notation and the results here, 
and discuss at qualitative level how the structure arises, whereas the details are discussed 
in Appendix |^. We shall also show concrete examples of the backgrounds as well as their 
/3- and 7-functions along the holographic RG flow below. 

The solid blue line represents a change in the UV asymptotics. The standard UV 
asymptotics of Appendix p.l.l| is obtained left of the blue line, whereas right of the blue 
line the solution "bounces back" at finite A, as discussed in Appendix D.1.2 . In the bounce 
back region the /3-function evaluated along the RG flow becomes zero at finite A, after 
which A starts growing toward the UV, and the standard UV singularity is not reached. 
Left the blue solid curve, the "standard" tachyon UV expansion of Appendix D.1.1 defines 
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the quark mass and the chiral condensate, which are not defined right of, or on, the blue 
fine. The red dashed fine has m = 0. This fine exists only for small x and terminates at a 
critical value Xc — 3.9959 (in scenario I), which matches with the definition of Eq. ( 8. IT] ). 
We stress that the solution with zero quark mass does not exist for x > Xc- Right of the 
red dashed line, the quark mass takes positive values and is monotonic in Tq. The contours 
give the quark mass, obtained by fitting the deep UV behavior of the tachyon solution to 
the expansion of Appendix The backgrounds in the contoured region will be identified 
as the physical ones, having lowest free energy. 

The value of the mass goes to zero when the solid blue (for x > Xc) or the red dashed 
(for X < Xc) curves are approached from within the contoured region as shown in Fig. |5|, 
which is not evident from Fig. ^ due to limited resolution. Between the blue solid and red 
dashed curves, the quark mass is small, but depends on x and Tq in a complicated manner, 
as we shall discuss below. 

The right hand plot in Fig. ^, which was obtained after modifying the potentials in the 
IR, shows similar qualitative features as the left hand plot. This is the case because the 
structure seen in Fig. |^ arises from the behavior of the solutions in the UV region and close 
to it, which is analytically tractable and almost independent of the change of the potentials, 
if, e.g., we keep the quark mass fixed (see the discussion in Appendix |f[) . However, the 
mapping to the IR asymptotics (in particular to Tq or ri ) is completely different in the two 
scenarios, which causes the differences between the plots. 

Let us then discuss how the structure of Fig. |^ arises from the background solutions 
(see Appendix ^ for a detailed analysis). The main point is that the closer we are to the 
thick blue curve (when approaching the curve from the right), the closer the background is 
to reaching the IR fixed point, when the tachyon finally grows large and drives the system 
away from it. Therefore the backgrounds near the blue curve will be quasiconformal, or 
"walking", so that the coupling A is approximately constant over a large range of ener- 
gies. The mass dependence can be then understood be studying the tachyon EoM and in 
particular the tachyon mass at the IR fixed point. Notice that as the red dashed curve of 
solutions with zero quark mass ends on the blue curve as x — )• Xc, quasiconformal behavior 
is expected in this region. 

To understand the behavior near the blue and red curves it is important to plot the 
UV parameter, m, versus the IR parameter Tq. We show this in figure |5| for x < Xc (left) 
and X > Xc (right). For x > Xc, there is a unique saddle point (regular classical solution) 
for each value of the quark mass. 

The situation for x < Xc is more complex and reflects the fact that the chiral condensate 
operator violates the BF bound in the (potential) IR fixed point. We see from the left of 
figure ^ that for each m > 0, there is a finite number of regular classical solutions that we 

^*The mass here is given in IR units, i.e., we actually plot mR = m/Aia, since we fixed the scale R of 
the IR expansions to unity in the numerics. 

^^As m and —m are related by a chiral rotation by tt, we expect that we can take m > 0. The chiral 
rotation is reflected in the background solutions in the symmetry T — >■ —T. Consequently, we can turn the 
negative mass solutions of Fig. |5|to solutions with positive mass \m\ by changing To — >■ — Tq. 

^"Notice that we have chosen potentials where the fixed point exist even at arbitrary small positive values 
of X. For another choice of potentials where the fixed point exist only up to a positive limiting value a;, (see 
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Figure 5: Left figure: Plot of tlie UV Mass parameter to, as a function of tlie IR Tq scale, for 
X < Xc- Right figure: Similar plot for x > Xc- The vertical solid blue and dashed red lines show 
where corresponding lines are intersected in Fig. ^. While these plots are for model functions, 
similar plots on how the bare quark mass m depends on Tq, for different x for the potentials of 
scenario I, can be found in figure |l^ in Appendix P. 



will label, by an integer n = 1, 2, . . ., with n = 1 being the rightmost solution in the figure 
(the one having the maximum Tq). 

Understanding of the qualitative shape of the left figure |5| is given by the fact that 
for X < Xc, T ^ 1 and A near the fixed point value A=k, the approximate solution of 
the tachyon is given by T ~ sin [A; log r + (j)]. Note that the constant k is fixed for 
fixed X but the normalization and (j) are determined by the boundary conditions and IR 
regularity. Therefore the tachyon starts at the boundary, evolves into the sinusoidal form 
for a while, and then at end diverges. Different solutions differ in the region in which they 
are sinusoidal, and it is this region that controls their number of zeros. This is explained 
in more detail in Appendix |^. 

For the solutions with high label n, tachyon changes sign several times before diverging 
in the IR. As we move to the left toward the vertical blue line in Fig. ^ (left), a new zero 
in the tachyon solution appears every time the mass curve crosses the horizontal axis. For 
m = we expect an infinite number of regular solutions for all positive integers n > 1. The 
presence of several such solutions reflect the violation of the BF bound, and are reflecting 
the Efimov minima seen in other contexts (see [0, |l^ ) . This agrees also with similar recent 
observations in fl^ . 

The hint of such multiple regular solutions / saddle points was seen already in [^] that 
treated the flavor sector in the quenched approximation. Indeed, a second regular solution 
was seen beyond the dominant one. In that case a calculation of the spectrum of mesons 
in this second solution indicated that this saddle point was unstable, as the spectra were 
tachyonic. 

It is interesting to point out that the presence of the Efimov-like tower of regular 
solutions is not tied uniquely to the existence of an IR fixed point solution in the landscape 
of the bulk theory that violates the BF bound. Even modifications of the potentials that 
do not allow this IR fixed point may still have the Efimov tower. The qualitative reason 

Sec. the structure is expected to be the one described here at least for a;, < x < Xc, and the structure 
in the region x < x, can be analyzed numerically. 
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Figure 6: The backgrounds with identicahy vanishing tachyon for x — 2 (left) and x = 4 (right). 
The red sohd, and bhie dashed curves are the vahies of A, and A as functions of log(r/^), respectively. 



is that once the bulk theory has a regime that is near critical, this is enough to trigger 



the presence of such multiple saddle points. This is indeed the case in |55| where in the 
quenched approximation a second solution exists even though there are no such IR fixed 
point. We have checked that exactly the same happens here in the probe limit a; — )■ 0, 
where the fixed point is absent, if the potentials of scenario I are used. Moreover, for the 
scenario II the full Efimov tower remains even in the probe limit. 

There is a more detailed analysis of the regular solutions in appendix ^. The compar- 
ison of the free energies of the various regular solutions is made in section 

8.5 Background solutions at vanishing quark mass 

To summarize the results of the above analysis, we identified the solutions of the contoured 
region of Fig. Q as the physical ones. For x < we found solutions for all m > 0, whereas 
for X > we found that m > so that the solution for m = was absent. 

We will discuss now how the background solutions vary as we move around Fig. Q 
(left). We will start with the massless case, where also an additional special background 
exist (for all x), for which the tachyon is identically zero. 

8.5.1 Solution with identically vanishing tachyon 

Let us start by analyzing the special solution with vanishing tachyon, which is not included 
in Fig. In the absence of the tachyon, the system is otherwise the one studied in |3f:], but 



with the dilaton potential Vg replaced by V^s = ^g~ xVfQ. This potential is guaranteed to 
have a maximum, corresponding to an IR fixed point, in the Banks-Zaks region since we 
matched it with the QCD /3-function. For the choice of Eqs. ( 8.24D -( |8.28| ) the fixed point 
actually exists for all < x < 11/2?^ 

There is a single solution to the equations of motion that reaches the IR fixed point, 
described in Sec. of Appendix ^. It is similar to the backgrounds studied in |44| 



where a /3-function inspired by supersymmetry was used. Similar backgrounds were also 
studied at finite temperature in [^]. We identify this special solution as the background 
corresponding to the chiral symmetry conserving phase, as the vev a vanishes. It is easy to 



^^This is also the case for the potential associated to scenario II, described in appendix]^. We could also 
slightly modify the potential so that the fixed point, disappears for x < Xm < Xc- We have also analyzed 
such a case and find only minor differences from those analyzed in detail in this paper. 
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Figure 7: The background for vanishing quark mass for various values of x (see the labels). The 
red solid, blue dashed, and magenta dotted curves are the values of A, A, and logT as functions of 
log(r/€), respectively. The thin lines are the UV and IR expansions of the solutions. 



construct the background numerically by shooting from the vicinity of the IR fixed point, 



e.g., by using the expansions of Appendix |E.3 . 

We plot the background for x = 2 and 4 in Fig. ^, where we fixed the scale R of the 
IR expansions ( [E.48| ) and ( E.4g| ) to one. The geometry is expected to asymptote to AdS 
both in the UV and in the IR, reflecting the flow from the IR fixed point to the standard 
UV one. Actually A is very closely linear in logr for all r, so that the deviation from 



AdS is not visible in the plots. Similar observation was made in |44|. While the solution 
with vanishing tachyon exist for all x, we will show in the next section that the other 
massless solution which involves a nontrivial tachyon profile and therefore chiral symmetry 
breaking (red dashed line in Fig. ^ has lower free energy whenever it exist, i.e., for x < Xc- 
Therefore, this background which correspond to a field theory flowing to an IR fixed point, 
is the physical one (for massless quarks) only for x > Xc, and Xc is indeed the edge of the 
conformal window. 



8.5.2 Solutions having nontrivial tachyon dependence 

The backgrounds having vanishing quark mass lie on the red dashed line in Fig. ^. We 
plot the corresponding background as a function of r for a few values of x in Fig. |^. As 
we matched with the IR expansions and chose their scale R to be unity, the IR scale 
is approximately fixed to 0{1). The changing of the background as the critical value 
Xc — 3.9959 is approached, is best visible in the solutions of A (the solid red curves). 

The dependence of the solutions on x meets the expectations from field theory. For 
X = 2 the solution is "running": A has simple and smooth dependence on logr. As x is 
increased to 3, a small distortion appears which becomes better visible for x = 3.5. The 
solution of A is developing a plateau at A ~ 25, as it approaches a fixed point. Indeed for 
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Figure 8: The /3-functions for vanishing quark mass for various values of x (see the labels). The 
red solid, blue dashed, and magenta dotted curves are the /3-functions corresponding to the full 
numerical solution (dX/dA) along the RG flow, the potential VcS — Vg — xV/o, and the potential 
Vg, respectively. 



X = 3.9 the coupling constant A "walks": it takes an approximately constant value as r 
changes by a few orders of magnitude. Such walking backgrounds have been studied in the 



context of IHQCD by using a model beta function in the probe limit in |46]. 

We also note that A (dashed blue curves) depends linearly on logr up to the IR 
region so that approximately Ac=l — log(r/£), and the metric is thus very close to the AdS 
one, even over the quasiconformal region where A walks. The tachyon (dotted magenta 
curves) is small and decoupled from the evolution of A and A in the UV and in the walking 
region. It becomes 0{1) (the curve crosses zero) only after the coupling has already 
started to diverge. This agrees with our expectation that the UV behavior, the behavior 
in the walking region, and in particular the phase structure is basically independent of the 
choices of IR behaviors of the potentials and the form of the tachyon action for large T. 

We also show the UV and IR expansions of the various fields, given in Eqs. ( p.2| ), ( p. 9] ) 
of Appendix^, and in Eqs. 8.29| , respectively, as thin lines where possible. These lines are 



often poorly visible since they overlap with the background. In the running region [x = 2) 
the full solution can be obtained to a good approximation by interpolating between the 
expansions. As x increases, the region of validity of the UV expansions is pushed to smaller 
logr, and neither UV nor IR expansions work in the walking region that grows as x — )• Xc. 
It is also illustrative to discuss the behavior of the system in terms of the /3-functions. 



^^It is difficult to obtain solutions with the value of mass close enough to zero to produce the (quadratically 
vanishing) tachyon dependence of the massless solutions when shooting from the IR due to limitations from 
numerical precision. Therefore, we matched the background solution at logr ~ —8 with a tachyon solution 
that was obtained by shooting from the UV and assuming decoupling, and plotted a combination of these 
to obtain the truly massless tachyon profile. A similar procedure was required for Fig. hd below. 
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Figure 9: The factor exp(— aT^) plotted against A along the RG flow for various values of x. The 
red solid, blue dashed, magenta dotted, and green dotdashed curves have x = 2, 3, 3.5, and 3.9, 
respectively. 



For this we recall that in the absence of the tachyon, potentials and beta functions are linked 
according to Eqs. ( ^.18 ), and (|8.19|) . By using this fact we can quantitatively estimate the 
validity of the tachyon decoupling in the UV and in the IR that was discussed above. First, 
recall that Vf{\,T) vanishes exponentially for large T. Therefore, in accordance with the 
discussion of Sec. |6|, the behavior of A and A is described by Vg{X) — xVfQ{\) {Vgil)) in the 
UV (IR) where T — )• (T — )• cxd). Eq. ( ^.19D gives directly the approximate beta function 
in the UV, whereas in the IR we must replace V^fr(A) — )• V^(A) as in the Yang-Mills case 

We show in Fig. ^ the /3-functions corresponding to the UV (dashed blue curves) and 
IR (dotted magenta curves) potentials, obtained by solving Eq. |8.1S| , and compare them to 
dX/dA evaluated along the RG flow of the numerical solution (red solid curves) for various 
values of x. First, notice that the x dependence of the effective /3-function dX/dA is as 
expected from Fig. |^: for x = 2 it is qualitatively similar to the Yang-Mills /3-function, 
and as x approaches Xc we find a typical quasiconformal behavior where the fixed point is 
almost reached at a finite value of the coupling. 

As A — )■ the effective /3-function dX/dA matches very well with the expectation from 
the tachyon decoupling (the dashed blue curves). Similarly, toward the IR (A — t- oo) the 
asymptotics of the red curves are similar to the magenta ones, which were obtained by 
taking A — t- oo. However the convergence towards the decoupling limit (blue curves in the 
UV, and magenta ones in the IR) is slower in the IR than in the UV. The main reason for 
this is understood by studying Fig. |^ Since we plot the /3-functions as functions of the 
coupling A, and A diverges in the IR much faster than the tachyon, the tachyon decoupling 
as A — 7- oo takes place slowly. This is confirmed in Fig. ^ by plotting the factor exp(— aT-^) 
of the tachyon potential, which controls the tachyon decoupling, as a function of A. Indeed 
this factor approaches the constant value of one quickly in the UV, whereas the convergence 
to zero in the IR is slower. 

Notice the clear similarity in the A dependencies of this factor and the /3-function 
dX/dA in Fig. |. 
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Figure 10: The gamma functions along tlie RG flow for vanishing quark mass as a function of 
A (left) and logr (right) for various values of x. The red solid, blue dashed, magenta dotted, and 
green dotdashed curves have a; = 2, 3, 3.5, and 3.9, respectively. 
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Figure 11: Left: The ratio of the UV (Auv) and IR (Air) scales as a function of x for m/Auv = 
10^^, lO"'^, . . . , 10 from thin to thick lines (top to bottom). The vertical dashed line is at critical 
Xc — 3.9959. Right: The ratio of the UV and IR scales as a function of m/Auv for x = 2, 3.5, 3.9, 
and 4.25 from thin to thick lines (bottom to top). 



Finally we plot the effective gamma function 

7 _ 1 dT _ dlogT 
T ~ TdA ~ dA 



.30) 



along the RG flow against A and logr in Fig. |To[ When the quark mass is zero, 7/T ap- 
proaches —3 in the deep UV (see Appendix D.l.lj ). In the UV region 7/T is approximately 



independent of x and increases with A until it reaches —2 near the value A ~ Ac where 
the fixed point develops as x — )• Xc- This is in line with the discussion of the preceding 
sections. In particular, when plotted as a function of logr we see that 7/T is close to —2 
in the walking region, corresponding to the saturation of the BF bound. We have checked 
that this behavior gets more pronounced as we choose values of x even closer to Xc so that 



a plateau near the value —2 develops in the left hand plot of Fig. |10 



8.6 Backgrounds at generic quark masses 

We now discuss the solutions of Fig. ^ for generic quark masses. As pointed out above, 
the massive solutions (the contoured region left of the blue and dashed red lines in Fig. ^ 
exist for all values of x. Except for the modified UV asymptotics of the tachyon, no new 
classes of qualitatively different backgrounds with respect to the massless case are found. 



-43- 



However, adding a mass introduces a new scale to the system, which affects the back- 
ground in a different way depending on its size and the value of x.^^ This effect may be 
illustrated by studying the ratio of the UV and IR scales of the background, which mea- 
sures how close to the IR fixed point the system comes. We define the scales in term of the 
UV and IR expansions, i.e., Auv = A in Eqs. ( p.2| ) and Air = 1/R in Eqs. ( ^.29 ), which 
we have fixed to one. 

The dependence of Auv/ Air on the quark mass in UV units m/Ayy and on x is 
depicted in Fig. |ll|, and meets the expectations from field theory. The left-hand plot shows 
the ratio as a function of x for various choices for the quark mass. We see that there is 
a qualitative difference between the regions with x < Xc and x > Xc- For x < Xc, chiral 
symmetry breaks spontaneously even for m = 0, and there is some range of small masses 
where the background is essentially independent of m. This is best seen on the right hand 
plot, where Auv /Air levels for small masses for values of x below the critical line (lowest 
curves). 

When the mass grows large enough (essentially larger than the scale of the spontaneous 
symmetry breaking), it starts to fix the IR scale directly, and Auv/Air decreases. For 
X > Xc there is no spontaneous chiral symmetry breaking, and the IR scale is determined 
smoothly by the value of the mass. Prom the log- log plot on the right we see that the 
dependence of Auv/Air on m is a power-law. Naively one could expect that the IR scale 
is directly given by the quark mass, m ~ Air, so that Auv/A-ir oc (m/Auv)~^ (which 
is also the result one gets by approximating T{r)/l = mr and using the arguments of 
Appendix |^). The nontrivial energy dependence of the quark mass modifies the power 
from —1 to smaller values. 

9. The free energy 

We now analyze the free energy for zero quark mass. In this case we identified two distinct 
solutions, one with identically vanishing tachyon and the other with nontrivial tachyon 
background. We shall show that the latter one is energetically favorable in the region 
where it exists (x < Xc). We start with the generic definition of free energy for our action. 
The free energy is given by the on-shell Euclidean action plus counterterms (which we 



will not need in this article). From (5.5) the Euclidean action takes the form: 



-(M^TV^) / d^x 



-X Vf{\ T) Vdet {gab + h{X, T)daT a^T) 
By using the Einstein equations ( |5.7| ) we can eliminate R, which leads to 

5os = -{M^N"^) j d^x 

X e>^ h{X, T) r'2 



(9.1) 



^ e'^ Vg{X) + ^ X e^^ Vj{X, T) ^l + e-2^MA,r)r'2 (9.2) 



3 Vi + e-2^/i(A,r) r'2 



'See ||69(| for an analysis within a different framework. 
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Next, one can solve Vg from (^^) and ( |5.9D : 



e^^ Vg{X) = i3A' e^^y + | e^^ Vf{X, T) I + /i(A, T) r'2 



+-e^^^^li£^=, (9.3) 



and inserting this into ( |9.3D the on-shell action can be integrated: 

5os = 2 iV2 y d^x e^^y = 2M^N^ V4 [A' e^^] [° . (9.4) 

Here e and rg are the UV and IR cutoffs, respectively. In all the backgrounds which we 

consider here the contribution from the IR vanishes, so we will drop that term. Following 
1 36] the Gibbons-Hawking boundary term is given by 

Sgh = 2M''N'^ [ d^xVhK =8M^ N'^V4A'i€)e^^^'\ (9.5) 

JdM 



where we have used that \3G\ K = 4e "^A' . We obtain the following expression for the free 
energy of the system: 

£ = 6M^N^V4 A'{e) e^^^'^ . (9.6) 
9.1 The free energy difference of the m = backgrounds 

The free energy calculated above is obviously divergent as e — )■ 0, and needs to be regu- 
larized. However, different solutions with the same UV boundary conditions (the quark 
mass m and the UV scale Auv) have the same divergent terms as well as counterterms and 
differ only through a finite term which can be extracted from the UV expansions. This is 
the case in particular for the two backgrounds having zero quark mass, one with vanishing 
tachyon and the other with a nontrivial tachyon solution, we now discuss in detail how the 
free energy difference between these two backgrounds can be obtained, by expanding all 
quantities as series at the UV singularity r = 0. Since the leading UV free energy behaves 
as corrections 0{r^) to the behavior of A and A will possibly contribute to the finite 
terms. 

As discussed above, in the UV the tachyon decouples from the equations of motion for 
A and A. For m = 0, the leading corrections to these equations due to the tachyon are 
suppressed by T"^ or e~2^T'2, i.e., by 0{r^). Therefore the coupling to tachyon does not 
affect the free energy directly in the massless case, and we can set it to zero. 

Let us take 

Kff(A) = ^^(A) - xVfiX, 0) = VgiX) - xVfoiX) = ^ [l + UiA + ^sA^ + • • • ] . (9.7) 

The finite contribution to the free energy can be studied by writing 

A{r) = Ao{r) + r^Ai{r) + 0{r^) (9.8) 
A(r) = Ao(r) + r%(r) + 0(r6) 
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in close analogy to Eqs. (^.S]), where Ai{r) and Aj(r) have now series expansions in 1/logr 
at r — 7- 0. The expansions of for and Aq are given in Eqs. ( p.2| ). Inserting these as well 
as suitable Ansatze for Ai and Ai in the EoMs of A and A, and expanding up to ©(r^) we 
find 



Ai{r) = A 



1 



19 



+ 



121og(rA) 

445^1^ - 320^2 
36V^ log(rA) 



(9.9) 



-5 + 



where A is the same UV scale that appears in the expansions of Aq and Aq, and A is 
a free parameter. It is recognized as an integration constant of the EoMs that did not 
appear in the leading order {0{r^)) UV expansions. Two solutions having the same A and 
m = 0, but different IR behavior, will have different A as its value is not fixed by the EoMs 
asymptotically in the UV. Inserting the results for Ai and Ai in the expression (9.6) for 
the free energy, we find 

AS = 6M^N^VifAA , (9.10) 

where AS and AA are the differences between the two solutions in free energy and the 
constant A, respectively. 

It is useful to also study the corresponding variation in the /3-function. We will actually 
use the phase variable defined by 



X(r) 



A'(r) 



3A(r)A(r) 



3A 



We may again write 



X{r) = Xo{r) + r^Xi(r) + 0{r^) . 
Substituting the expansions ( |9.8D , (|9.9D in the definition of X, we identify 

Xi{r) = A 



(9.11) 



(9.12) 



15, , 5[80yi2_64V2-log[-log(rA)] (23^12 -64^2)] 



■ log(rA) + 



+ 



(9.13) 

By using the logarithmic expansion from above, the result can be expressed in terms of A: 

"log(9Vi/8) (23^12 _ 64V2)" 



on 4 

X(A) = Xo(A) + -=^exp 



9^2 



32 14_64V2 



3^1 A4 

Finally, we recall that when the tachyon has decoupled, X satisfies 



aa:'(a) 



+ 



3A d 
X{X)dX 



logoff (A) 



X{Xf 



1 



6 



[1 + 0(A)] (9.14) 



(9.15) 



Indeed it is easy to check that the result ( |9.14 ) is consistent with this equation, and that 
A is the integration constant which parametrizes all solutions of the differential equation 
for given potential V. 

The remaining task is to extract the coefficients A from the chiral symmetry breaking 
(with nontrivial tachyon profile) and conserving backgrounds (with identically vanishing 
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|A£|/A^, 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



Figure 12: The free energy difference between the chiral symmetry conserving (T = 0) and 
breaking (T — oo as r — > oo) sohitions. Sohd blue curve is the difference for the chiral symmetry 
breaking solution having a monotonic tachyon with no zeroes, and the dashed red curve is the 
difference for the solution having one zero. 



tachyon) which have zero quark mass, and then use the formula ( |9.10D to calculate the free 
energy difference. Extraction of the coefficients is done by studying the variation Xi of the 
phase function, and details are given in Appendix p.l . 



We find that the chiral symmetry breaking solution is the energetically favorable one 
in the region where it exists, i.e., for < x < Xc- We plot the free energy difference (setting 
M'^N'^V4, = l) in Fig. ^ (soUd blue curve). 

Notice that approaches zero both as x — )• and In the case x — )• 0, 

we do expect that the effect vanishes as indeed the number of flavors, which controls the 
backreaction of the tachyon, vanishes. One expects linear dependence oc x: since the 
background configurations behave smoothly as x — )• 0, the energy difference arises due to 
the explicit x dependence in the action and due to its linear effect on the background. The 
case X ^ Xc will be discussed in the next section. For Xc < x < 5.5 only one solution with 
zero quark mass exists (the one with T = 0), as discussed in Sec. so there is no need for 
comparison. 

There are also solutions with zero quark mass where the tachyon has one or more 
zeroes (see Appendix ^). We have verified that these solutions have larger free energies 
than the one without a tachyon zero. The red dashed line in Fig. 12 is the free energy 
difference between the background with vanishing tachyon and the one with a nontrivial 
tachyon solution having one zero. The energy difference between the T = solution and 
the solutions having more than one tachyon zeroes is even smaller. 

Therefore, for x < Xc and m = 0, the standard tachyon solution has the lowest free 
energy and the chirally symmetric one, T = 0, has the largest free energy. All other 
undulating solutions have free energies that are between these two (and closest to the 
T = solution). 

For \m\ > 0, the standard tachyon solution has the lowest free energy and the non- 
standard ones, have higher free energy. 



-47- 



10. Scaling below the conformal window 



10.1 BKT scaling of the chiral condensate 

We shall now argue that the chiral condensate oc a obeys the BKT 
scaling behavior, 

a oc exp I -= 

\ yJXc — 

as X approaches the critical value Xc where the solution ceases to exist. This argument will 



or Miransky B 



(10.1) 



be supported with numerical results in section 10.3. The ratio of the UV and the IR energy 
scales will show similar scaling. This behavior is known to arise both in Dyson-Schwinger 
fo], ^] and holographic approaches [0, 11, 46 1. Indeed, our analysis has many similarities 
with both Dyson-Schwinger and earlier holographic approaches, and in particular with the 



recent study [12| of a related model. 

We shall not give a precise proof but only sketch how the scaling arises. It is enough 
to study the region near the UV where the tachyon is small r(r) <C 1, so the tachyon 
decouples from the EoMs of A and A. We will neglect the logarithmic corrections to the 
tachyon (see Appendix p. 1.1] ) which play no role in the scaling argument. 

In the deep UV, where the coupling is small A ^ 1, the tachyon behaves as 



T(r) 



or 



(10.2) 



As r increases T stays small, and A starts to approach the fixed point value A = A,,, which 
maximizes V^(A) — xV^^iX). The behavior of A and A in this region is given by the T = 
asymptotics of Sec. 



A = A, + O 
A = - log(r - ro) + + C 



ruv 



-2S 



(10.3) 
(10.4) 



where 5 is a positive parameter defined in Eq. ( [E.51| ) . This approximation is valid for some 
intermediate region ruv ^ ^ '^iRj where at the scale hr the tachyon finally becomes 
0(1), and drives the system away from the IR fixed point. We choose x near the critical 
value, < Xc — X ^ 1. The tachyon IR mass was calculated above in Sec. 



where 



Air (4 - A 



G(A,x) 



ir) = -miViR = G(A*,x) 



24a(A) 



(10.5) 



(10.6) 



/i(A)(yg(A)-xl7o(A)) 
Keeping formally x fixed while varying A*, the right hand side will become equal to 
four as A* reaches a critical value Ac: G(Ac,x) = 4. The scaling behavior will appear as 
A* —5- Ac from above, and then also x — >• Xc- We expand around this point: 

d 



Air(4-Air) = G(A,,x) 



4 + ^G(Ac,x)(A, 



Ac) + --- = 4 + k(A*-Ac) + 



(10.7) 
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For < A* — Ac ^ 1, we obtain 

Air ~ 2 ± - Ac) . (10.8) 

The tachyon solution becomes 

T{r)^Qp(—] sin (^/k{X, - AJ log — + A . (10.9) 
When A moves away from the fixed point, while moving towards the UV, it will at some 



point become smaller than Ac so that the asymptotic solution ( |10.9| ) fails. However, as we 
shall see, the leading scaling originates from the region where A,,, — A ^ A=k — Ac, and we need 
the solution for smaller A only to make contact with the deep UV behavior and the definition 
of a. Even for A < Ac, k(Ac — A) is small for almost the whole range of r with ruv < ^• 
Therefore, within the this region, it is sufficient to neglect k(Ac — A) and approximate 
T(r) oc r^. Further, we introduce an intermediate scale f where A=k — A ~ A* — Ac, and 
write the approximation as 

T{r)^c(—) ; ruv«r<f, A<Ac (10.10) 

T{r) ~ Qp ( — ) sin ( ^'^(A* - Ac) log t + ; r < r < hr , Ac < A < A, (10.11) 

where we use r instead of ruv ^ the reference value of the logarithm for later convenience. 

For r > riR there is no obvious way to write a good approximation for the tachyon 
solution. However, as we shall see, such an approximation is not necessary for finding the 
scaling behavior. 

The scaling behavior can be found by matching the tachyon solutions in the different 
regions. First, we require that the solutions of Eq. (10.2) and ( 10.10| ) join approximately 
continuously at r ~ ruv- This gives 

C^ar^v- (10.12) 

Notice that a is not a free parameter here, but it will later be fixed by the matching 
procedure. Further requiring approximate continuity at r ~ f we find^^ 

C7fp~ar^v ; (10.13) 

The remaining task is to match with the unknown IR behavior at r ~ rjR. First, hr was 
defined as the scale where the tachyon becomes 0{1) and drives the system away from the 
fixed point, so r(riR) ~ 1, which fixes riR in terms of a: 

Qp ~ <Trl^ ~ f — ) • (10.14) 

Finally, we need to match the solution to T'(riR), which is ©(I/hr), since the tachyon 
EoM is apparently regular in this region. Notice that we must indeed fix this number 



*This matching involves a subtlety which does not affect the scaling (see Appendix I 
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to have a solution that asymptotes to the "good" singularity in the IR. The good IR 
asymptotics have one free parameter, the normalization of the tachyon in the IR. This is 
however aheady determined by requiring T{riji) ~ 1. Therefore, the argument of the sine 



function in (10.11) is basicaUy fixed to a given 0(1) number at r = hr, which gives the 
desired BKT scaling: 

V'^(A* - Ae) log ^ = 0(1) (10.15) 

so that 

with K positive. Finally, we notice that the connection between f and ruv can be obtained 
from ( 10. 3| ) by using the definition of f. This results in a power law: 



(A, -Ac)"^ (10.17) 

ruv 

which can be neglected as a subleading correction to the exponential scaling. Taking this 
into account 

The inverse of this scaling result is expected to hold for any ratio of IR and UV energy 
scales independently of their precise definitions. By using Eq. ( |10.14 ) we find the scaling 
results for a: 

1 / 2K \ 1 f 2k \ 

^uv V \/A* — Ac y r^v \ — X j 

Notice that Xc and Ac were defined by G(A=i, (xc), Xc) = 4 and G(Ac,x) = 4, respectively, so 
that A* = Ac at X = Xc- Since G is smooth in this region. A* — Ac could readily be replaced 
by Xc — X in the results above (after rescaling K). 

In the above discussion we neglected several subtleties. These issues are analyzed in 
Appendix |H| In particular, after more careful analysis, we are able to find explicit results 
for K and K: 

K = ^= . ; K = = . (10.20) 

^ ^JmGiXcx) ^-|:G(A.(x),x)|_^ 

As a final remark, we stress that the above sketch was to a large extent independent of 
the details of the model. In particular, we did not need any information on the nonlinear 
terms in the tachyon EoM and on how the IR boundary conditions are fixed. 

10.2 Scaling of the free energy 

Let us then study the scaling of the free energy difference A£ between the solutions without 
and with chiral symmetry breaking (at zero quarks mass), which was studied numerically 
in Sec. |9|. Linearizing the EoMs for A and A by writing 

A{r) = Ao{r) + r^Ai{r) + 0{r^) (10.21) 
A(r) = Ao(r) + r%(r) + 0(r6) , 
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r 



Figure 13: The tachyon logT (left) and the coupHng A (right) as functions of logr for an extreme 
walking background with x = 3.992. The thin lines on the left hand plot are the approximations 
used to derive the BKT scaling (see the text for explanation), and the vertical dashed (dotted) lines 
mark logruv (logr). 



we related A£ with the variation of the leading coefficients of Ai and Ai when expressed 
as a series in logr at r = 0. The source of this variation is the difference in the tachyon 
solution, and it can be analyzed in the limit x — )• Xc- The tachyon contributes corrections 
©(T^) and {e'^^{T'f) to the EoMs of A and A. For zero quark mass these contributions 



are 0{r^) in the deep UV and thus decoupled from the variations of Eqs. (f0.21). However, 
as we have pointed out, in the walking region ruv ^ ^ ^ Hr the tachyon is ©(r^), and 
therefore the tachyon contributions ©(r^) do couple to Ai and Ai.^^ Consequently, we 
expect that, e.g., 

r^AAi{r)^T{rf ; ruv < r < hr , (10.22) 

where AAi{r) denotes the difference in Ai between the solutions with intact and broken 
chiral symmetry. We expect that to a good enough approximation, the size of Ai equals 
the coefficient A defined in Eq. (^^) even in the walking region where the UV series does 
not converge. Therefore, we obtain 

Ai ~ ctV^v (10.23) 
and finally the scaling result for the free energy difference reads 

AS .,00 , / AK \ , I AK 



Ai~ASv~r„-*exp(-^^)~.,;texp , (10.24) 



The next-to-leading terms of A and A are still ©(r") in this region since there are no sources which 
could change this. 

^^This can be verified by inserting the Ansatze (f0.2f ) to the EoMs and studying the solutions similarly 



as done for Xi in Appendix S.l 
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Figure 14: Left: log(cr/A'^) as a function of x (dots), compared to a BKT scaling fit (solid line). 
The vertical dotted line lies at a; = Xc- Right: the same curve on log-log scale, using Ax ~ Xc — x. 
To guide the eye we added the straight dashed line corresponding to the BKT scaling fit without a 
constant term, that was fixed to go through the data point with smallest Ax. 

log(Auv/AiR) 




Auv/AiR 



Figure 15: Left: log(Auv/AiR) as a function of x (dots), compared to a BKT scaling fit (solid 
line). Right: a / plotted against Auv/Air on log- log scale. The line is a fit with the expected 
power law. 



10.3 Comparison with numerical results 



We now compare the analytic results above to numerical solutions. In Fig. 13 we plot the 
tachyon (left) and the coupling A (right) as functions of log r in an extreme walking case 
with X = 3.992 such that Xc — x ~ 0.004. The numerical tachyon solution with zero quark 
mass (blue thick curve on the left) was obtained by gluing together the various solutions 



described in Appendix G.2, We compare the solution to the analytic approximations of 
Eqs. (10.2), ( 10.10| ) and ( |10.11 ), shown as thin green dotted, magenta dashed, and solid 
red curves, respectively. The parameters of these curves were chosen such that a has the 
extracted value (see Appendix G.2| ), ruv = with Auv obtained by fitting A to its 

UV expansion, and f is the value where A reaches Ac. The parameters Ac, A,,, and k were 
calculated directly from the potentials, and (j) was given an arbitrary small value. The 
agreement between the approximation and the full numerical solution is remarkably good. 
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We have also compared the expected scahng to the values of chiral condensate that 
were extracted from the background (see Sec. ^ as detailed in Appendix p. 2 . The results 
for various values of x are the dots in Fig. ||. The soHd Hne is a fit to the BKT scaling 
behavior, given by 

a 6.8 



8.6 



(10.25) 



which describes the data well. The analytic result of Eq. ( 10.20D gives for the potentials 
used 

2k = 6.10 (10.26) 

which agrees with the fitted result 6.8 within the precision of the fit. We have checked that 
using only a few of the data points with the highest x brings the fitted valued of K closer 
to the analytical one. 

In Fig. 15 (left) we compare the ratio of the scales Auv = ^ and Air = 1 / R as defined 
by the UV and IR expansions of the background, respectively, to a BKT scaling fit with 
K ~ 3.4. We also checked that the scaling of Eq. ( 10.14 ), 



A3 



njv 



VAuvy 



(10.27) 



is satisfied to a high precision on Fig. |l5| (right). 



11. Conclusions 



We analyzed a novel class of holographic models (V-QCD), which reproduces the main 
features of QCD in the Veneziano limit of large Nf and Nc with x = Nf/Nc fixed. 

V-QCD is on one hand based on a successful holographic model of Yang-Mills (YM) 
theory, and termed improved holographic QCD (IHQCD). IHQCD contains a dilaton cou- 
pled to five-dimensional gravity background. Its characteristic feature is the holographic 
renormalization group flow of the YM coupling constant, identified as the exponential of 
the dilaton, as a function of the energy scale, identified roughly as the inverse of the bulk 
coordinate. On the other hand, the model builds on earlier work on including matter in 
holographic models via flavor branes in the quenched approximation, i.e., neglecting the 
backreaction of the brane on the dilaton and the background metric. In particular, we use 
the tachyon Dirac-Born-Infeld (DBI) action originally introduced by Sen. 

Putting together these two frameworks, the dynamics in the Veneziano limit is modeled 
by a system of a dilaton and a tachyon coupled to five-dimensional gravity. The dilaton 
action is fixed as in IHQCD. For the tachyon we consider a generalized DBI action, where 
the dilaton dependence is parametrized in terms of a few potentials, which are a priori 
unknown. We may then constrain the unknown potentials, among other methods, by re- 
quiring that the UV physics implements that of (perturbative) QCD and that the solutions 
are regular in the IR. 

The essential observation for uncovering the dynamics of the system, is the identifica- 
tion of the effective potential. It involves terms from both the dilaton and DBI actions, and 
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takes the role of the dilaton potential of IHQCD. For large x, the perturbative Bank-Zaks 
IR fixed point can, and must be implemented trough the effective potential. Even further 
away from the Banks-Zaks region, the solutions can continue to flow to the fixed point in 
the IR. Interestingly, the fixed point can also be "screened" by the tachyon dynamics, such 
that the theory comes very close to it and the coupling almost freezes, but eventually starts 
running again in the deep IR. This kind of backgrounds are termed as the quasiconformal 
or "walking" ones. 

In this article we did a detailed analysis of the backgrounds and the zero-temperature 
phase structure of V-QCD. Our main results are as follows: 

• We generalize the holographic RG flow of IHQCD to include the evolutions of both the 
dilaton and the tachyon, which are controlled by the holographic /3- and 7-functions, 
respectively. Remarkably, for potentials that are analytic in the UV, the interplay of 
the dilaton and the tachyon automatically results in the anomalous dimension of the 
quark mass having physically reasonable UV asymptotics, i.e., a power series in the 
't Hooft coupling. 

• If we require that the UV expansions of the potentials capture the essential QCD 
physics, and choose potentials that join smoothly with their UV expansions, the V- 
QCD phase diagram is always the physically relevant one. That is, for zero quark mass 
we find a phase transition at a; = Xc- The conformal window, where the backgrounds 
have an IR fixed point, extends from a; = Xc up to the maximal value x = 11/2 where 
asymptotic freedom is lost. Below x = Xc, chiral symmetry is broken, and the theory 
has similar behavior in the deep IR as the QCD we have observed in Nature. 

• The edge of the conformal window is stabilized such that the dimension of the quark 
mass at the IR fixed point approaches two (the anomalous dimension approaches 

one) as the edge is approached. Under reasonable assumptions for the potentials, we 
find values of Xc within a narrow band around the value Xc = 4. High value of the 
anomalous dimension is of importance for applications to (walking) technicolor. 

• Below Xc but close to the edge of the conformal window, we find quasi-conformal or 
walking backgrounds. 

• Backgrounds for (any) nonzero quark mass exist. In the conformal window the quark 
mass triggers chiral symmetry breaking, and below Xc introducing the quark mass 
affects the scales of the theory in a physically reasonable manner. 

• We verify that the standard vacua which give the above physical phase structure 
also have the lowest free energies. They have monotonically increasing or identically 
vanishing tachyon profiles, when the chiral symmetry is broken or intact, respectively. 
For X < Xc we also found a tower of unstable "Efimov" vacua, for which the tachyon 
can change its sign arbitrary many times. 

• Finally, as x approaches Xc from below for the backgrounds with zero quark mass, 
we show that the chiral condensate, as well as the length of energy range where the 
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background stays close to the fixed point, obey the characteristic Miransky or BKT 
scahng law. 
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APPENDIX 



A. General solutions for /3 and 7 

It is not difficult to construct numerically j3 and 7 which solve the partial differential 
equations (|7.5|) and (7^). We can use the fact that they were derived from a set of 



ordinary differential equations. More precisely, the equation have the form 

/3^+7^ =-^/3(/3,7) (A.l) 

/3g+7^ =-^7(/3,7) (A.2) 

where and J%y are independent of the partial derivatives of the /3- and 7-functions (and 
the dependence on potentials was left implicit). 

Now we can apply a standard method for solving first order partial differential equa- 
tions (PDEs), which in this case applies even for a pair of PDEs. We first look for curves 
along which the system reduces to ordinary DEs. Not surprisingly, such curves coincide 
with the holographic RG flow. That is, we can implicitly define such family of curves, by 
requiring that each curve (A(A),T(A)), where A parametrizes the curve, satisfies 

A'(^) = I3{\{A),T{A)) ; T'{A) = j{X{A),T{A)) . (A.3) 

Then we notice that the differential operators in ( |A.l| ) and (|A.2| ) become derivatives along 
the curves: 

^f^{X{A),T{A))=Tp{P{A)MA)) ; ^j{X{A),T{A)) = T,{(3{A)MA)) ■ 

(A.4) 

Further, the equations along the curves are essentially in one-to-one correspondence with 
the system (5.8)-(5.11). Indeed, the original system of DEs can be formally recovered be 



eliminating the /3 and 7 functions by using Eqs. ( [A.3| ). Therefore, any solution of the 
original system satisfies the PDEs along the curve (A(r),r(r)) that it defines. 

As the final step of the method, we notice that the PDEs only depended on the 
derivatives of /3 and 7 along the curves. The derivatives perpendicular to the curves 
can be freely chosen. Therefore, any (continuously parametrized) family of curves which 
satisfies ( [A.3| ) will define a solution to the PDEs in some region of the (A, T)-plane. That 
is, the general solution to the PDEs is given by the planes that the solutions of Eqs. (| 



( 5.11 ) draw in (/3, A,T) and (7, A, T)-spaces as the boundary conditions to the equations 
are varied in an arbitrary manner. 

Remarkably, the amount of degrees of freedom of the solution matches with the expec- 
tation for a system of two first order PDEs, for which the solution should depend on two 
arbitrary functions. Since /?, 7, A, and T are all invariant under the symmetry of (p.4|), 
the number of integration constants in the system ( ^.8D -( pT| ) which are relevant for the 
solution of PDEs is three. 

A generic one-dimensional family of these parameters, and consequently a generic 
solution to the PDEs, can be defined by giving the dependence of two of these in terms of 
the third one which makes two arbitrary functions. 
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We have demonstrated how the boundary conditions of the PDEs are mapped to those 
of the original set of EoMs. Now the analysis of the solutions of the EoMs (see Appendices ^ 
and ^) suggests at least two natural ways to choose special one-parameter families of the 
boundary conditions, and to define special solutions to the PDEs. 



1. Require that the solutions to EoMs end in the "good" IR singularity of Sec. [E.2.2 
(with varying quark mass). 



2. Require the more generic IR singularity of Sec. E.2.1 and keep the quark mass fixed 



(In the second case, the "good" IR singularity is expected to arise as some limit of the 
more generic IR behavior so that it appears at a boundary of the region in (A, T)-space 
where the P and 7 functions are defined.) 

B. A coordinate transformation 

It turns out that it is convenient to solve the system numerically using ^ as a coordinate 
instead of r. Therefore, we present the system after this transformation: 



3A2 



l + g-2/i(A,T) r2 



A , 4A q\ A2 ?, ^ I dVg I -2 dVg , xt^ dh 



A A g A A2 8 \ dX V 2g2 Jl+q-^h ' 



f + At--t 



A + X -^logiVhVf) 



5r °^ 



-TAAiog(.y,)4^iogy, 

where the dots are derivatives with respect to A and we defined 

,(^) = e-|^. (B.2) 

C. Examples of explicit potential choices 

We will construct explicit examples of potentials, which give physically reasonable back- 
grounds. We consider an Ansatz for Vf of the form 

y/(A,r) = v>o(A)e-"W^' , (C.l) 

and assume that h does not depend on T. First, we expect in the UV 

Vg{X) = Vo + ViX + V2X'^ + --- (C.2) 
VfiX,T = 0) = VfoiX) = Wo + WiX + W2X'^ + ■■■ (C.3) 
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The potentials must also produce the good kind of IR singularity discussed in Appendix E.2.2 
with P = 1/2 |]3^ , which constrains the asymptotics of Vg to 



Vg{\) ~ A^/=^7b^ ; A 



(C.4) 



(assuming that V/o plays no role in the IR for the chiral symmetry breaking solution as 
tachyon and/or A diverge so that exp(— a(A)T^) tends to zero). 

In addition, the IR behavior of V/o should be chosen such that the fixed point (max- 
imum) of V^(A) — xV/o(A), which is guaranteed to be present at large x — )• 11/2 if we fix 
the small A series of the potential appropriately, continues to exists up to sufficiently low 
X. This is most easily achieved if Vyo diverges faster than as A — )• oo, so that the fixed 
point exists for all x. Another possibility, mentioned in Sec. is that the fixed point exist 
only for x > x*, where x=k > is relatively small. We have checked that choosing such 
potentials does not change any results at qualitative level, and do not discuss this choice 
further here. We also require that the potentials are analytic at A = 0. A simple Ansatz 
that meets these requirements and involves free coefficients up to two-loop order is 



Vfo{\) 



Fo + ^lA + V2 



A^ 



(l + A/Ao)2/3 



Vl + log(l + A/Ao) 



Notice that the AdS radius 



12 



(C.5) 
(C.6) 

(C.7) 

must be well defined for all x up to x = 11/2, which sets an upper bound for Wq for given 
Vq. We also expect Vf to be positive at small A.^"^ 
Therefore, we take 

< 1^0 < . (C.8) 
However, as discussed in Appendix p.l.l| , the sum of the anomalous dimensions of the quark 



mass and the chiral condensate is not equal to 4 if Wq = (then 5 = 1 in Appendix D.1.1). 
Therefore, we discard this option. Notice also that if we saturate the upper limit with 
Wq = ^Vo, the AdS radius diverges in the Banks-Zaks limit x — )• 11/2 unless we choose 
an x-dependent Vq. 

In addition, we can fix the coefficients Vj, Wi by mapping to the field theory /3-functions 

as 



Vg{\) - xy/o(A) 



12 
12 



8,YM. , 23(6^)^-366™ 



l + -6o-n + 
1 + -60A + ° 



81 



366i 2 



81 



A^ + • • • 



(C.9) 
(C.IO) 



Here the (scheme independent) QCD /3-function in the Veneziano limit with vanishing 
quark masses up to two-loop order are 



^^Negative Wo is also problematic since, at least for small |Wo|, it generates a zero of Vf, which causes 
the tachyon solution to become singular. 
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2 11 - 2x 



2 34 - 13x 

'3 (4^)4 



and = bi\x=o- Setting Vq = 12 and Wq = 12/11, which hes in the middle of the allowed 
region of Eq. ( p. 81 ) , we obtain 

44 _ 4619 

9^ ' 

4(33 - 2x) 



Vi -- 

Wi 



V2 
W2 



38887r4 ' 

23473 - 2726x + 92x2 



(C.12) 
(C.13) 



997r2 ' 427687r4 

Further, we choose Aq = 871^ to prevent the higher order terms in the UV expansion of the 
potentials from growing unnaturally large. 

In addition, we need to choose the functions a(A) and h{X) appropriately. As discussed 
in Section |8| and in Appendix p.l.l| , we must have 

h{X) 2f 



a{X) 



(C.14) 



Here the coefficient hi can be matched with the one-loop anomalous dimension of the quark 
mass, which reads in the Veneziano limit 



7m (A) 



(47r)- 



rA + 



By matching with Eq. ( D.l[l| ) from Appendix, we obtain 



(4vr) 



from which 



48, 4 ■ 

On H /ll 

39 3 



115 - 16x 



2167r2 



(C.15) 



(C.16) 



(C.17) 



The IR behavior of the functions h{X) and a(A) is linked to the tachyon behavior in 
the IR for the solution which breaks chiral symmetry. They are discussed in detail in 
Appendix E.2.2| , where essentially only two different cases, which are consistent with the 
good IR singularity, are found. These possibilities are produced by the following choices. 



I We can choose the function a(A) to be constant, and the function h{X) to have power- 
law IR asymptotics: 

1 3 _ 3 



h{X) 



(1 - ^A)4/3 



a(A) = ao 



22 



;ii-x) , 



(C.18) 



which corresponds to the special case of p = 4/3 and = in Appendix E.2.2. In 
this case the tachyon diverges exponentially. 



r(r) ~ To exp 



81 3^/6(115 - 16x)^/3(ii -x) r 
812944 2V6 R 



(C.19) 



as r — )• 00. Here R is the IR scale of the solutions which is defined by Eqs. (EJi 
( [E.20| ) and (E.21). Tq is the only free parameter. 
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II The other choice is practicahy a generalization of the first one. It has a simple form 
of h{X), but more complicated a(A): 



h(X) = a(A) = -(ll-x)V-^^±^ , (C.20) 

^ ' (l + A/Ao)4/3 ' ^ ^ 22^ ' (i + A/Ao)^/3 ' v ; 

so that p = 4/3 and a = 2/3 in Appendix [E.2.2| . The extra term proportional to /12 
was added to make a positive. We choose its coefficient to be small, /12 = l/Ag- The 
tachyon behaves as 

^ ^ \ R ^ ' 

for large r. Here ri is a free parameter. 

We have checked that both scenarios lead to qualitatively similar results. In the nu- 
merical calculations we use, for definiteness, the choice I, unless stated otherwise. 



D. UV behavior 

In this Appendix we shall discuss the UV behavior of the system ( ^.8D -( |5T1 ) in general. 



This analysis should be compared to that carried out in |3G, ^] in the absence of the tachyon 
backreaction. Recall that apart from the two degrees of freedom of the transformation (^^) , 
the solutions contain three integration constants. In the discussion below, the degrees of 
freedom refer to these three "nontrivial" constants. 

We shall not discuss the most general behavior of the solutions, but make some phys- 
ically motivated assumptions. In particular, we restrict ourselves to the potentials Vg and 
Vf which are bounded as A — )• 0, and which are smooth at any finite A. In general we 
are interested two types of potentials: ones that start from a constant value at A = 0, are 
monotonic as A increases, and approach -|-oo as A — )• 00, and ones that start at a constant 
value at A = 0, increase until they reach a maximum at some A = A,,,, and thereafter 
monotonically decrease to —00 as A — )■ cxd. This should be kept in mind while reading 
the analysis below, as some of the arguments below may fail for more generic potentials, 
even though no assumptions are listed explicitly. Further, we mostly restrict to effective 
;9-functions /3efr = dX/dA which are negative. 

D.l Generic behavior 
D.1.1 Singularity at A = 

Near the standard UV singularity, the tachyon must behave schematically as T(r) ~ mr + 
ar^ whereas the other fields have a logarithmic dependence on r. Therefore the tachyon 
decouples asymptotically as r — t- 0, and we may analyze its UV behavior by first solving A 
and A with T = and then by analyzing the tachyon EoM for this background. 
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Asymptotic behavior of A and A We take 



12 



^eff(A) = Vg{X) - xVfiX,0) = — [l + Vi\ + V2X^ + 

Then the (leading) UV expansions of A and A can be written as 

4 



A{r) 



log - + 



1 

162 



91og(rA) 

+ m log [- log(rA)] 



(D.l) 



(D.2) 



95-^ 



-23+^ 



ViA(r) 



log(rA)2 
log [- log(rA)] 



91og(rA) 



+ 



46 _ 128 V2 
81 sTT^ 



log[rA]^ 



+ 



+ 
1 



log(rA)3 



log(rA)'^ 



Notice that they contain no free parameters (in addition to A). In fact, after using the 
equations of motion there is one degree of freedom left in the coefficients of the above 
expansion, but as it turns out, this freedom can be eliminated by rescaling A. We have 
removed this extra parameter by requiring that the coefficient of the l/(logrA)^ term in 
the expansion of A vanishes. 



Tachyon UV asymptotics We take 



Vf{X,T) = e 



-a{\)T'^ 



Vf{\) 



and parametrize 



12 



^eff(A) = Vg{\) - xVf{\ 0) = [1 + ViX + VsA^ + 
xVf{\) = A^ \Wq + WiX + W2\^ + ■■■] 



h{\) 2e' 



a(A) 



— \l + hiX + h2X^ + ■■■' 
6 



(D.3) 



(D.4) 



(D.5) 



where 5 is a nonnegative integer. Here the leading coefficient of h/a was already fixed in 
order to have the correct UV mass of the tachyon . We further assume that 



(D.6) 



/i(A) = A«(1 + 0(A)) . 
It is enough to study the linear terms in the tachyon EoM, which read 



T"{r) + 
3 + 



+ 



log(rA) 
8{hi+Vi) 



+ 



log(rA)2 



+ 



T(r) 



T'{r) 



. 



(D.7) 



3Filog(rA) Vlog(rA)2 
From this one could expect that the solution has the form 

T(r) ~ mr (1 + O (1/logr)) + ar^ (1 + O (1/logr)) 



(D. 
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i.e., the logarithmically suppressed corrections to the EoM show up as logarithmically 
suppressed corrections to the functions. However, this is not the case: an Ansatz for the 
solution which assumes this kind of corrections fails. The correct asymptotics reads 



1 4 s i , iH 

-r(r) = mr(-log(rA))» 2 2 + 3^1 

+ar^{-log{rA)) '^^^ 2 + 2 ' 



^ ^ C7i + C72log(-log(rA)) ^ 1 



log(rA) \log(rA)2 
^ ^ I)i + Zj2log(-log(r-A)) ^^f 1 



(D. 



log(rA) \log(rA)2 



where Ci and Di are known functions of 6, ^, hi, Vi, V2, Wi, and W2, which we suppressed 
for brevity. 

The "surprising" logarithmic power corrections in Eq. ( p. 9] ) can actually be identified 
as the nontrivial running of the quark mass and the condensate in the UV, which arises 
as their anomalous dimensions are different from zero. To make this explicit, we calculate 
the gamma function T' /A' in the UV. For m 7^ it is dominated by the linear tachyon 
solution: 

T TA' log(rA) \\og{rAy J ^ ' 

whereas for m = we find 

4 _ 3(3+g) , 4hi / 1 \ 

J = _3 + ^_^_±J^ + 0(_i_) . (D.ll) 



T log(rA) Vlog(^^)^ 



The next-to- leading terms in ( |D.1C| ) and ( p. 11 ) are mapped to the one-loop anomalous 



dimensions of the quark mass and the chiral condensate in QCD, respectively. Since they 
should add up to zero, we must have 5 + =0. The easiest way to satisfy this is to take 
(5 = = ^. In particular, the expression h{X) = A~^/^ with ^ = —4/3, which was found in 
the probe limit |3^, does not work, since 5 was required to be an integer to ensure that 
the /3-functions have power series with integer powers at A = 0. 

Finally, it is easy to verify that the UV expansion presented here match with those of 
Sec. which were derived by using the holographic beta functions. 

D.1.2 A bounce back at finite A 

A bounce back may take place when the potential Uofr(A) = V^(A) — xV/(A,T = 0) has a 
maximum at some A = A,, signaling the presence of an infrared fixed point, and U'(A) < 
for A > A=K. If the tachyon is sufficiently small, the effective /3- function dX/dA hits zero, and 
becomes positive when the system is evolved toward the UV. Therefore the coupling has a 
finite minimum (> A*) and the above "standard" UV singularity at A = is not reached. 
All fields are analytic at the point where A' = 0. The bounce back behavior is found in the 
white regions of Fig. Examples of the /3 and 7-functions evaluated along the RG flow 
for the bounce back scenario are shown as the dotted curves in Fig. 16 in Appendix ^. 



D.2 Special case: UV fixed point at finite A 

We have identified one special UV singularity, which is found as a limiting case between 
the two first generic behaviors discussed above (the blue curve of Fig. In this case 
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the asymptotic solution is expected to depend on two integration constants. The solution 
terminates as the /3- function dX/dA approaches zero at the maximum A* of the effective 
potential Ves{X) = Vg{X) — xVf{X,T = 0). The singularity is found at a fixed value of r = r* 
where A diverges and A approaches A* from above. Examples of the /3 and 7-functions 
evaluated along the RG flow with this UV fixed point are shown as the thick blue curves 
of the middle and right columns of Fig. 16 in Appendix ^. 

Depending on the value of the tachyon mass at the fixed point (see Sec. 8^), the 
asymptotics may be written in two different forms. Let us recall the definition of the 
dimension A at the fixed point: 

24a(A,) 



A(4- A) 



(D.12) 



h{X,)V,s{X,) ■ 

When X > Xc the right hand side of the definition is smaller than 4 so that there are 
two real roots A = A±. The geometry approaches the AdS one near the fixed point, 

A{r) = - logr + log 4 + Ao{Arf^- + ■■■ (D.13) 
A(r) = A* + Ao(Ar)2^- + ••• (D.14) 
r(r) = To(Ar)^- +••• , (D.15) 

where A_ is the smaller root, = 12/14ff(A*), we chose r* = 0, and the constants Aq, 
Aq, and To satisfy two constraints which can be solved from the EoMs (5.8)-(5.11). There 
are two free parameters which can be taken to be the coefficients of the tachyon solutions 
with the dimensions A±. The solution associated to A_|_ will appear at the next-to-leading 
order only if we choose the coefficient Tq of the above solution to vanish. 

When X < Xc we have two complex roots A± = 2 zizik. Now the asymptotics reads 



A{r) 

A(r) 
T(r) 



log r + log 4 + (Ar-)'' Ai + Ai sin {2k log(Ar) + 0^) 



+ ■■■ 



X^ + {Ar)^ Ai + Ai sin (2A: log(Ar) , 
ro(Ar)^sin(A;log(Ar) 0t) H ■ 



+ ■■■ 



(D.16) 

(D.17) 
(D.18) 



The coefficients in the next-to-leading terms for A and A can be solved from the EoMs 
( |5.8| )-( pT ) by inserting the tachyon asymptotics. The free parameters are Tq and (pT in 
this case. 



E. IR behavior 



In this Appendix we discuss the IR behavior of the system ( |5.8D -( ^T1 ) in general. As above, 
we restrict to certain quite simple potentials and to cases where the /3-function dX/dA takes 
negative values in the vicinity of the IR singularity. The IR structure is much richer than 
the UV one mostly because there are much less obvious constraints on the potentials. 

We shall discuss here both generic and special singularities. Of these the most in- 
teresting ones for us will be the the special ones that depend only on one free parameter 
(excluding trivial reparametrization symmetries). They can be identified as the "good", 
or fully repulsive, IR singularities. Examples of such singularities are identified below in 



Sec. E.2.3 partially based on the analysis of Sec. E.2.2. 
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E.l Generic cases 

There are two generic IR "singularities" which in fact do not involve divergences of any of 
the fields. Therefore, they also appear independently of the details of the potentials to a 
large extent. 

E.1.1 Divergence of the derivative of the tachyon 



A typical, generic IR behavior is similar to what was found in the probe limit in |36|, where 
the tachyon goes to a constant value but its derivative diverges. Indeed, the Ansatz 

T{r) = n+Ti./¥7^ + T2in-r) + --- (E.l) 
A{r) = Aiin - r) + A2{n - rf^^ + ■■■ (E.2) 
A(r) = A, + Ai(r, - r) + A2(r, - rf^^ + ■■■ (E.3) 

solves the equations of motion ( |5.8D -( |5?Tl| ) quite in general. We do not present the rather 
complicated constraint equations which follow for the constants in the expansions, but it is 
not difficult to check that the solution has three independent integration constants and is 
thus indeed generic. Since none of the fields diverge at r = r^,, it is natural to take all the 
potentials to be analytic at the point of expansion, and the solution is expected to exist to 
a large extent independently of the choices for them. 

Notice also that there is no real singularity at r = r^,: one can make a coordinate 
transformation such that all fields are analytic in the vicinity of this point. A natural 
choice that realizes this is to use T as the coordinate. Of course, this leads to all fields 
being double- valued functions of r, with the two branches having the same absolute value 
of Ti but opposite signs. If this is allowed, it is not hard to find analytic solutions which, 
for example, start at a UV singularity at A = 0, bounce back at a point where the tachyon 
derivative diverges, and return to another singularity at A = 0. 

E.1.2 A bounce back as d\/dA — > 

There is also a three-dimensional space of solutions where the coupling reaches a maximum 
value, and then starts to decrease with decreasing A so that the /3-function d\/dA becomes 
positive. All fields are analytic in r at the point where the /3-function is zero. We have not 
checked how the solutions continue to evolve in the region of positive /3-function. 

E.2 Special IR singularities 

In addition to the generic IR behaviors discussed above, we have identified several true 
IR singularities, where the fields A and A diverge. Note that above, the divergence of the 
tachyon derivative was proportional to the parameter Ti, which could take both positive 
and negative values. Solutions with arbitrary small |Ti| also exist, as well as the limit 
Ti — 7- where the solution, which in general ends in the divergence of A and A rather than 
T. Therefore we expect that spaces of generic solutions with positive and negative Ti, 
respectively, will be separated by a subspace of solutions with singularity in the IR. The 
behavior of the system in this case depends strongly on the asymptotics of the potentials. 
However, the tachyon often decouples asymptotically from A and A, in particular if the 
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tachyon diverges in the IR, which is the expected behavior for physicahy relevant singu- 
larities 55]. Here we shall assume that the decoupling takes place, since completely 
general classification of the singularities seems daunting. 

As the tachyon decouples, the classification of singularities for A and A follows earlier 
studies |3^, We shall review the results here for clarity. Assuming that the tachyon 
tends to Tq in the IR, the effective potential that drives the metric and the coupling in the 
IR is 

Vm{X) = Vg{X)-xVf{X,To) (E.4) 
where Tq can be infinite in which case T^r(A) = V^(A). We parametrize 

y,a(A) = (log A)- (y„ + + ^ + . . . ) (E,5) 



as A —7- oo. The equations of motion are Eqs. ( |5.8| ) and (5.9) with T(r) = Tq so that 
T'(r) = 0. There are two types of singularities (see p8[): "generic" ones where 



in the IR, and "special" ones where X — )• — |Q. 
E.2.1 Generic metric singularity 

The generic singularities exist for Q < 4/3 and depend on one free parameter. The system 
is solved by the Ansatz 

A= - log 5r- log R + Ai6A^-^'^^/^{- log 6r)'' + ■■■ (E.7) 
3 

A = -L(^Ao + Ai5r(8-'5«)/3(-log5r)^ + ---) (E.8) 

where 

6r = {n-r)/R (E.9) 

is the "conformally invariant" distance from the singularity, Aq is the free parameter, and 
the dropped terms are suppressed by 1 / log 6r. Plugging this in the equations of motion, 

_ 27(2Q - mxl^+' 
^ 16(44- 57Q + 18Q2) " ^^"'^^ 

For Q = 4/3 the singularity exists if P < 1. In this case the asymptotics reads 

A = - log 5r -logR+ ° _^ (- log 6rf+^ + • • • (E.12) 
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E.2.2 Special metric singularity 

This kind of singularity exist for < Q < 4/3 so that the asymptotic value —4(5/3 of 
X lies between zero and one. We need to require that the potential is asymptotically 
positive, Vo > 0. The solution does not involve any integration constants in addition to the 
ones linked to the reparametrization symmetry, which suggest that when combined with a 
proper tachyon solution, "good" IR asymptotics can be identified. 
If2/3<(5<4/3we find a singularity at finite value of r: 

1 



A 



9Q74- 1 

Q 



log 5r — \ogR + O (1/log 5r) 



(E.14) 



g^-4/9^°^'^ + ^ 



2 log 2 + (1 - 2P) log 3 + (P - 2) log (9Q2 _ 4) (E.15) 



+ log (16 - 9Q2) - P log Q - P log(- log 5r) - log Vq 



+ 0{l/\og5r) 



where again 6r = (r^ — r)/R. Here R and are the integration constants which reflect 
the reparametrization symmetry, but no other free parameters appear. 
If < Q < 2/3 similar formulas hold for r — )• 00: 



A 
X 



1 



1 - 9QV4 

Q 



logr — log R-\-0 (1/log r) 



(E.16) 



4/9 -Q2 2Q 



2 log 2 + (1 - 2P) log 3 + (P - 2) log (4 - 9^^) (E.17) 



log(l6-9Q2) - PlogQ - Plog(logr) - logFo 



0(1/ logr) , 



where now r = (r — r^j/R. 

If Q = 2/3, and P < 1, there is a singularity at r = 00. The asymptotic solution reads 



A 



r -rp 
R 



IP, R 5 P 13 2^1 

+ ^0 - i:-. 5 log + 7: + T + o^log- + 

21 — F r — ro 6 4 2 2 oPVq 



(E.18) 



+- 

log A = + 

+ 
where 



-52P^V^ + 4P^V^ + 27P^V^ + - 64PV^ + I28PF0V2 / R 



288P(1 + P)V^ 



r — ro 



3 f r — tq 



5 3P 
4 ~ ~8 



2 \ R J 4 8 PVq 

2opVo2 - mp^v;^ + 9P^yo^ - uv^ + upv^ - i28pyo^2 / R 



(E.19) 



192P(1 + p)v;} 



r -ro 



+ 



1 



a 



R 



1 - P 

2^3^-^ 
(l-P)e^o^ • 



(E.20) 
(E.21) 



As pointed out in [^], this special case produces a good match with the IR physics 
of QCD, in particular if we choose P = 1/2. We will confirm below that the potentials 
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of the tachyon action can be chosen such that the tachyon asymptotics also meets all 
requirements known to us, and the produced singularity is of the "good" kind. We will 
use these singularities in our analysis, and the above formulas will be used to fix the IR 
boundary conditions for the numerical solutions. 

If Q = 2/3, and P > 1, we find a singularity at finite value r = r,,, of the coordinate. 
The asymptotics transforms to 

+^0-9l slog + ? + T+9^1°g9+^ ^-22 

— r J 21 — F — r 6 4 2 2 SPVq 

-52PVo^ + AP^V^ + 27P^V§ + QWl - UPV^ + 128^^01^2 / r, -r \° 
^ 288P(1 + P)Vi \ R ) ^ " 

, , 3 / P \° 5 3P Vi ,^ , 

-20PVo^ - mP^V^ + QP^Fq^ - 64^1^ + UPVl - 128PVoy2 / r, - . 
^ 192P(1 + P)V^ V R 



where a = 1/(P — 1) and Aq is related to R as in Eq. ( E.21 ). 



Finally, for (5 = 2/3 and P = 1 the metric factor A diverges exponentially as r — )• oo, 

A = -exp(^)-logP + ^ + llog6-ilogyo + ^ + ^ (E.24) 
, 128^2-211/0 / r-ro\ , 

3 /r-ro\ 13 , -51^0-128^2 / ^"^oV ..^ 

E.2.3 Tachyon behavior 

To complete the analysis, one should insert each of the above asymptotics to the tachyon 
EoM and check what the tachyon asymptotics is for various choices of the potentials Vf 
and /i, and start looking for the "good" kind of singularities. Once the potentials are 
fixed, one can check if a solution, which is consistent with the assumption that the tachyon 
decouples, indeed exists. If it does, it can depend on one or two additional parameters. 
For the good, fully repulsive singularities the number of free parameters (excluding those 
related to the reparametrization symmetry) is equal to one, i.e., we must have a special 
metric singularity combined with a one-parameter tachyon asymptotics. In addition we 
should require that the tachyon diverges in the IR, since that kind of solutions have bulk 
flavor anomalies similar to those of QCD |^, 

Here we shall restrict to the special metric IR singularity with Q = 2/3 and P < 1, 
since it is expected to include the most interesting cases due to additional constraints from 
confinement and excitation spectra p6|. We parametrize 



V/(A, T) = VfoiX) exp {-a{\)T'^) ; h = h{\) (E.26) 
With this parametrization, the tachyon EoM reads 

T" + PiT' + P2T + PaT'^ + P4T'2t + F^T'T^ + F^^T'^T^ = (E.27) 
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where 



Fi = 3A' + A'— log(/i(A)y/o(A)) ; 



F3 = /i(A)e-2^ 



4yl' + A'— log(xA(A)y/o(A)) 



F2 
Fa 

Fa 



2a(A)e^-^ 
h{X) 

2a(A) ; 



(E.28) 
(E.29) 



and the primes are derivatives with respect to r. Notice that the last two terms vanish if 
a(A) is constant. 

We consider generic power-law asymptotics 

h{X) ~ hoX-P ; a(A) ~ aoA'" ; Vfo{X) ~ WqA^ (E.31) 

of the potentials at large A, and introduce a shorthand notation for the asymptotic behavior 
in pl8D -( F:20D : 



r -rp 
R 



logA = ^f^ 
^ 2\ R 



+ log + Ac + 



Xc + O 



R 

r — ro 
R 



r -rp 
R 



where we set rp = 0. Then the leading behavior of the coefficients Fi at large r is 
3a(2 -\- p — t) / r 



Fi 
F2 
F3 
Fa 
F5 



2r 

2aoe2^'=+^^'^+'^-^'= 
hp 



a-l 



/ioa(16 — 6r + 3p)e 



^) exp 



'2Ac-pXc 



3 3 

2^ + r-2 



2aoe''^'= exp 

Sapaae'^^'^ 
2r 

3ao/ioQ^'7e 



4i? 
3cr / r \ " 

7^ \ a 



^) exp 



-2Ac-pAc+o-Ac 



■ exp 



3 / r \ a 



R 
R 



2R 



exp 



3 3 
' + r-2^ 



(E.32) 
(E.33) 

(E.34) 
(E.35) 
(E.36) 
(E.37) 
(E.38) 
(E.39) 



For (7 = the leading terms of F^ and Fq, given above, become zero. If a(A) is 
constant for all A, these terms actually vanish identically. If a (A) only asymptotes to a 
constant value, the leading behavior of F^ and Fq is determined by the next-to-leading 
term in the asymptotics of a(A). These will contribute in some particular cases, as we 
discuss below. 

• o" > and p < 4/3. The tachyon EoM ( [E.27 ) is dominated by the terms oc T, T'T'^, 
which have the coefficients F2 and -F5, respectively. Solving the tachyon from these 
terms gives the asymptotics 



r(r) ~ To 



SaaTp 



dr exp 



-p-2 



(E.40) 
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Substituting this back to the full equation of motion, we see that with the above 
constraints for a and p it is indeed a solution. Further, the factor exp (— a(A)r^) 
vanishes double-exponentially, which confirms the decoupling of the tachyon from 
the other fields. In addition to the trivial reparametrization symmetry, the only free 
parameter of the asymptotic solution is Tq, which suggest that the singularity is of 
the "good" kind. However, tachyon solutions that are regular in the IR have bulk 
flavor anomalies which differ from those of QCD 55]. Therefore, we discard this 
option. 

cr > and p = 4/3. The same terms continue to dominate, but the asymptotic 
changes. We find instead 



Tir) ^2J^^^e^^+i^^^/F^i (E.41) 
V 3acrno 

where ri is a free parameter. One can again check that this is indeed a solution, and 
that the tachyon decouples from A and A. The terms proportional to T'^T and T'^T"^ 
are suppressed only by r~°, but taking them into account results in a trivial factor 
multiplying the equation of motion, so that the solution in Eq. ( |E.41 ) is unchanged 



for any value of a. The asymptotics has only one free parameter, and the tachyon 
diverges as r — )• oo, so this solution is acceptable. 

(7 < and p < 4/3 — a. We find two different cases. The asymptotics is qualitatively 
similar to ( E.40D , but arises in a slightly different way. The leading terms are those 



proportional to T' and T, corresponding to coefficients Fi and F2, respectively, as 
well as the double derivative term T" . This term being leading, one might expect 
that the asymptotics contains two free parameters, and is thus not of the good kind. 
This is indeed the case for r > 2 + a. For r < 2 + a only one parameter family of the 
asymptotic solutions is consistent with the other terms being subleading. However, 
since the tachyon becomes constant in the IR, we have the aforementioned problem 
with flavor anomalies, and hence we shall discard this solution in any case. 

o" = 0, p = 4/3, and r < 10/3. The leading terms are typically those proportional 
to T'^ and T''^T. However, if o(A) is not constant, there is an extra constraint 
from the next-to-leading term in the expansion of a(A). If, for example, a(A) = 
Oo -|- ai/X + • • • , we need to require a > 1. Assuming that all constraints are met, we 
find the exponential behavior 

r(r) ~ Toe^'' (E.42) 



where Tq is a free parameter and 



4 Ae+4Ac/3 n 

(10 — STjariQ 



This special case is the asymptotics that was discussed in the probe limit in |]36|, and 
is acceptable. By inserting the expressions for Ac and A^ from Eqs. (^1|), (p3o| ). 
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and ( [E.21 ), the coefficient simplffies to 



(10-3t)(1 -P)/ioFo« 

£7 = 0,/? = 4/3, and r > 10/3. The leading terms are proportional to T' and T, which 
results in the tachyon vanishing asymptotically. In this case the factor exp (— o(A)r^) 
goes to one rather zero, which suggest that the correct physical picture, as discussed 
in the main text, cannot be achieved, even though the tachyon apparently decouples 
in the IR. 

cj < and p = 4/3 — cr. The leading terms are again proportional to T' and T. For 
r < 10/3 — cj/2, the solution is exponentially increasing, 

r(r) ~ e^^ (E.45) 

with 

^ = (10 - 3r - 36)aho ' ^^'^^^ 

The factor exp (— a(A)T^) vanishes in the IR limit if a < 1, which is the expected 
behavior, so the solution is acceptable. For a > 1 or r > 10/3 — a/2 the factor goes 
to one instead, and it seems that the correct physical picture cannot be obtained. 
The borderline case a = 1 has either behavior depending on the values of other 
parameters. 



In all the remaining cases, in particular for large p, the asymptotic solution of Eq. ( E.27 ) 
oscillates with increasing frequency as r — )• oo. Therefore, the tachyon is apparently not 
decoupled from the other fields, and the singularity discussed in the subsection does not 
exist. 

In summary, we found good solutions only if p = 4/3 and cr > 0, or p = 4/3 — cr and 
cr < 0. The latter case is included for < a < 1, and the former case we found some extra 
constraints at the endpoint cr = 0, which are detailed above. In the numerics we shall fix 
P = 1/2 so that a = 2, and use potentials with p = 4/3 and a > 0, see Appendix ^ 

E.3 Singularity at the IR fixed point with T = 

As discussed in the main text, some of the solutions with identically vanishing tachyon are 
expected to correspond to field theories where the chiral symmetry is conserved. Potentials 



for Yang-Mills theory were discussed extensively in 35]. Therefore, we will only consider 
the case of potentials Veff(A) = V^(A) — xVf{X,0) which have a (single) maximum at some 
A = A*, interpreted as an infrared fixed point. 

In the absence of the tachyon the space of solutions is one-dimensional. We identify 
two distinct one-parameter families of solutions: one where the solution bounces back to 
smaller couplings as the /3-function dX/dA goes to zero (a special case of the bounce-back 
solutions discussed above in Sec. [E.1.2D and another where the /3-function asymptotes as 



dX/dA ~ — 3A (a special case of the solutions discussed in Sec. E.2.1| ). These families are 
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separated by a single solution where the /3-function terminates at zero at the maximum of 
the potential A = A,,, . 

The limiting solution has a singularity at r = oo where A diverges and A approaches 
A^, from below. We expand the potential around A = A* as 

Ves{\) = Vo + V2{X-Kf + --- (E.47) 

where V2 is negative. The equations of motion are solved by 

A = A.-(^) ' + ••• (E.48) 

A = -log(r-ro) + ^o + ^i f^-^j +••• (E.49) 

where ^0 and 5 are related to the IR AdS radius and the derivative of the /3-function, 
respectively, by 



4 = e^^» = ^ (E.50) 



and we also find 



U^_L_4|:>=^=, 4-2M-2, (E.51) 



''^ (E.52) 



9(25 -l)Ar 
E.3.1 Generalization to T <C 1 

A generalization of the T = case, which is of high interest to us, is where the tachyon 
mass and the condensate are very small in the UV, such that the tachyon remains small 
(|T| <C 1) even as A approaches the fixed point (A,,, — A <C 1). However, for any nonzero 
tachyon profile, the tachyon will eventually become 0(1) at some high r = rip, and drive 
the fiow away from the fixed point. Setting tq = 0, the above formulas ( E.4^ ) and ( E.4g| ) 



then hold in the limit R <^ r <^ rjR (with the understanding that depending on the value of 
S, the next-to-leading terms may be affected by the tachyon solution). This approximation 
is useful in the quasiconformal, or "walking" region of backgrounds for x below, but close 
to the edge of the conformal window at Xc- Notice that R is the scale where the coupling 
starts to deviate significantly from its fixed point value A^,, and may be therefore identified 
as the UV scale of the theory. Therefore we shall denote R = ruv in this case. 

The tachyon profile can also be derived in this region. Inserting the solutions of 
Eqs. ( p^ ) and ( |K49D into the tachyon EoM dsTn] ), and taking |r| < 1, we get 



T"{r) + -T'{r) + 7 7 f T(r) = . (E.53) 
r n(A*jr^ 

There are two kinds of solutions depending on the value of the tachyon mass term. As in 



Sec. p. 3 we can define the dimension A 

c 

7(A*) h{K)V,s{K) 



A(4-A) = ^"()-)f = 24a(A) 
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If the combination on the right hand side is less than 4, which corresponds to a; > Xc we 
find two real roots A±. In this case the tachyon solution is 

T(r) ~ To (—) (E.55) 

for rjR ^ r ^ rm , where A_ is the smaller root (unless we tuned the boundary conditions 
such that the quark mass is very small, in which case the solution with A_ — t- A-|_ needs 
to be included). 

For X < Xc we have two complex roots A± = 2 zizik instead, and the tachyon behaves 

as 

r(r)~To( — ] sin I fclog — + </. ) . (E.56) 

This oscillating solution is the root of the rich structure of backgrounds found for x < Xc 
in Sec. |8|, which are also discussed below in Appendix ^. 



F. Structure of the background as a function of x and Tq 

In this Appendix we will explain in detail how the phase structure of the background 
solutions seen in Fig. ^, and in particular the region with nearly conformal behavior, arises. 

The structure in the plots is linked to the transition of the system from the UV region, 
where the tachyon is small and the background is characterized by the potential T4fj(A) = 
V^(A) — xVfQ{X), to the IR region, where the tachyon is large, and the background is 
characterized by Vg{X). First, recall that Vef[{X) has a maximum at some A = A* which 
depends on x. For x — )■ we find from Eqs. ( |8.24[ ) and ( |8.26 ) that A* — )• oo, whereas for 



x —7- 11/2 we obtain A* — t- 0. This maximum suggests a presence of an IR fixed point of 
the /3-function for the coupling A. 

However, for a nontrivial tachyon profile the fixed point is not reached. When ap- 
proaching A = A=K from the UV, the solution is driven away from the fixed point as soon as 
the tachyon becomes large, T ~ 0{1), and the system enters the region where all EoMs are 
nontrivially coupled. The point where this happens, is controlled by the normalization of 
the tachyon in the IR, i.e., the value of Tq (Assuming potentials of scenario I). The larger 
Tq, the smaller is the value of A where the tachyon decouples. 

The solid blue curve in Fig. is the critical value of Tq where the tachyon decouples 



at A = A=K. Actually, precisely at this curve the UV asymptotics is that of Sec. |D.2| , and 
the solution ends at the fixed point. If Tq is smaller than the critical value, the tachyon 
will become small <C 1 during the flow toward the UV when we still have A > A*, so that 
the beta function corresponding to the effective potential Ves is positive. As the tachyon 
decouples, the beta function flow approaches that deflned by V^fj. Therefore, to the the 
region with "standard" UV behavior is not reached at all, but the solution bounces back 
at a flnite value of the coupling (see also Fig. |l^ below) where the beta function (evaluated 
along the RG flow) crosses zero. The strong dependence of the blue curve on x is explained 
by the dependence of A,,, on x: for example as x — >• 0, the fixed point moves to large values of 
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A*, and the critical tachyon IR value Tq required for avoiding the bounce back approaches 
zero. 

The argument above is not rigorous, as it involves the location of the decoupling of the 
tachyon which is not defined precisely. However, one should notice that as the blue curve is 
approached from above, the system is on the verge of reaching an IR fixed point so that the 
coupling freezes, i.e., it evolves very slowly for a large range of r. Meanwhile, the tachyon 
grows relatively fast with r. Therefore, the value of A where the tachyon decoupling takes 
place becomes more and more precisely defined as the blue curve is approached from above. 
This can also be seen in the numerical examples below and in Section |8[ 

To understand how the red dashed curve arises, we need to study the tachyon solutions 
in the UV region. Below the red curve the tachyon solution develops a zero, as is required 
by the negative value of the quark mass. This zero actually appears in the region, where 
(the absolute value of) the tachyon is still small. If we continue on the plot toward lower 
values of Tq, the quark mass becomes again zero, and then positive extremely close to 
the solid blue curve. This happens as the tachyon develops a second zero in the UV. We 
can continue further, and find solutions with an arbitrary number n of zeroes (which are 
very hard to construct numerically). See Fig. || (left) in the main text for the qualitative 
behavior of the quark mass as the blue curve is approached. 

The oscillating behavior of the tachyon in the UV is linked to the violation of the BF 
bound: when Air(4 — Air) is smaller than 4, the solutions for Air are complex, which 
results in the oscillations of the tachyon solution. We plotted the squared mass of the 
tachyon at the IR fixed point in Fig. ^, where the solid thick blue curve corresponds to the 
present choice of parameters. 

As we approach the solid blue curve from within the contoured region in Fig. 0, the 



background can be approximated near the fixed point as discussed above in Sec. E.3.1 



By the definition of Eq. ( |8.17D , the BF bound is violated at the fixed point for x < Xc- 
Therefore, as we approach the blue curve from above in this region and the system is 
about to develop a fixed point, the tachyon necessarily oscillates as soon as values of A 
close enough to A* are reached. In this case, the tachyon is well approximated by the 
solution of Eq. ( p3^ ), 

/ r \ " . / . . r 



r(r)~To — sin fclog— + , (F.l) 



where k is fixed by the potentials. 



^""f^-' i , (F.2) 



/i(A.)v;,(A.; 



but (p and Tq are free parameters, which will be fixed by the boundary conditions. For 
X > Xc, the BF bound is not violated at the fixed point, and therefore no oscillations 
are expected, and the tachyon dependence is as in Eq. ( E.55D . The (uppermost) curve of 



zero quark mass (the red dashed one in the plots) essentially limits the region of oscillat- 
ing tachyon: for example on the left hand plot above this curve quark mass is positive 
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(no tachyon zeroes) and below it the mass is negative (one tachyon zero), see the mass 
dependence in Fig. ^. Notice that the red curve must therefore join the blue curve at 

As we approach the fixed point keeping x, the range ruv r <^ rjR where Eq. ( |F.1| ) 
holds grows without limit. Here ruv is the scale where the coupling starts to deviate 
significantly from its fixed point value A=k as we follow the flow toward the UV. According 



to Eq. (F.l), the number of encountered tachyon zeroes is 



vrn ~ k log . (F.3) 

ruv 

We may take one step further, and find the scaling of the quark mass and the chiral 
condensate as ni^/ruv ~^ oo and n — )• oo. This can be done by matching the "intermediate" 
tachyon solution (F^l) with the UV (and IR) solutions (see Sec. 10 and Appendix ^ where 



we do the matching procedure more carefully for x — >• Xc)- The tachyon is supposed to 
become large at r ~ rjR, so Tq ~ 1. On the other hand at r ~ ruv we enter the standard 
UV region, where roughly T ~ mr + ar^. Matching in the UV gives the typical sizes for 
m and a for large n: 

3 / ruY\'^ f 27rn\ 
mruv ~ o-'^uv ~ i 7^ ) ^ i 'k' ) ■ ^ 

In particular, the maximal masses for which solutions with n tachyon nodes exist, or in 
other words the sizes of the bumps in Fig. |5| (left), numbered from right to left, must obey 
this scaling law. As the mass scale vanishes exponentially for n — )• oo, for any fixed m ^ 
there are only finite number of backgrounds as n is limited from above, whereas for m = 
we find an infinite tower of solutions. 

At this point is good to remind that the solution with no tachyon nodes (n = 0) always 
has smaller the free energy than the solutions with n > (see Sec. |9|). Also, the n = 
solution is not found in the scaling region where the system is close to having a fixed point 
in general. For m = this solution (red dashed curve of Fig. ^ enters the scaling region 
(which is close to the solid blue curve in the same figure) only in the limit x — )• Xc, which 
will be discussed in detail in Sec. 

We conclude with one more observation on the location of the curve of the vanishing 
quark mass. The above discussion was relying on the dimension Air at the fixed point, and 
therefore we could argue how the tachyon behaves in the limit where the background comes 
arbitrarily close to the fixed point, i.e., as we approach the blue solid curve from above in 



^*We cannot prove analytically that the scaling region is accessible, because this would require a de- 
scription of the solution in the IR region where tachyon is not decoupled. However, we may solve the 
EoMs starting from the UV with arbitrary small quark masses and vevs, which are guaranteed to enter the 
oscillating region, and have arbitrary many tachyon zeroes. The tricky issue is, if we can tune the solutions 
such that it ends in the good IR singularity after the tachyon finally grows large and the oscillations end 
in the IR. Due to the oscillating nature of the solution, it is plausible that the good IR singularity is found 
after any number of oscillations n, if it is found, e.g., for small n. In the end this question is settled by 
the numerics, which supports our expectations, and solutions with regular IR behavior are indeed found 
for any n. 
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Fig. ^. However, we can also present rough qualitative arguments on the behavior of the 
curve further away from the solid blue curve. Because the tachyon is in any case small in 
the interesting region, we can read directly from the linearized tachyon EoM whether it 
oscillates or not. The EoM reads 



T"(r) + 



3A' + X'^log{h{X)Vfo{X)) 



2A 



T'{r) + =^ • 



Recall that the tachyon is decoupled in this region and evolves independently of the 
other fields. Assuming decoupling the coefficients are essentially independent of Tq and 
can be solved directly from the potential Vc^[\) = Vg{\) — xVfQ{\). A is quite well 
approximated by ^ ~ — logr,^^ the second term in the square brackets of the coefficient 
of T' is small, and therefore the essential term is, as in the case of the fixed point, the 
ratio a(A)//i(A), which increases with A. As A grows, at some critical Ac the ratio becomes 
large enough, and the tachyon starts to oscillate. In terms of the gamma function, this 
means that 7/T reaches the value of approximately —2. The oscillations do not take place 



if tachyon grows large already for A < Ac so that nonlinear terms contribute in Eq. F.5 
Therefore, for the limiting solutions, the tachyon becomes 0{\) roughly at A = Ac- As the 
start of the oscillations means that the quark mass goes to zero, this mechanism also fixes 
the location of the curve with zero quark mass (red dashed curve in Fig. ^ . 

In summary, the solid blue curve of Fig. ^ is stabilized by the tachyon growing large 
at A = A*, whereas the dashed red curve is stabilized by the tachyon growing large at 
A = Ac- The two curves meet when Ac = A=k, which gives an alternative way to formulate 
the definition of Xc- 

It is interesting to compare the picture in our model to that arising in the Dyson- 
Schwinger approach in the rainbow approximation (see, e.g., [^]). Also in this framework 
it is useful to define two values of the coupling, corresponding to A=k and Ac above. The 
definitions are similar as here: A=k is the zero of the /3-function, and Ac is the value of A 



where the anomalous dimension of the chiral condensate reaches unity (so that A of Sec. 8.3 
equals two). 

Indeed the latter definition corresponds to saturating the BE bound in the present 
approach, which in the vicinity of the IR fixed point means the start of the tachyon oscil- 
lations, matching with our definition of Ac. Similarly as in our model, Ac = A* at the edge 
of the conformal window. A similar description of the conformal transition was also found 
in the holographic model of 

F.l Numerical analysis 

We illustrate the analysis above by studying the background numerically near the red 
dashed and blue solid curves of Fig. g 

We choose as reference values x = 2, 3.9 and 4.25, which have qualitatively different 
behavior. For zero mass, they correspond to field theories with running, walking and IR 



conformal behavior of the coupling constant. Fig. 16 shows, for the above values of x, the 



^^Explicit expressions can be derived in some approximation schemes, e.g., as series in the limit of small 
/3-function. 
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Figure 16: The dependence of the background on Tq near the red and bhie hnes of Fig. |[ Top 
row: X = 2. Middle row: x = 3.9. Bottom row: x = 4.25. Left column: the quark mass as a 
function of Tq. The vertical blue solid and dashed red lines mark the solutions terminating at an 
UV fixed point, and having m — 0, respectively, as in Fig. ^j. Middle column: /3-functions dX/dA as 
Tq is varied over the range of the mass plots (left) with constant steps, shown as thin black curves. 
Dotted, dashed and solid curves are /3-functions for backgrounds with a bounce back towards the 
IR, m < 0, and to > 0, respectively. The limiting cases between these behaviors are the given by 
the thick blue solid curves, which terminate at an UV fixed point, and the thick red dashed ones, 
which have m = 0. The magenta dotted curve is the /3- function when the tachyon is completely 
decoupled, solved from Eq. ( ^.19 ). Right column: gamma functions 1/T x dT/dA as To is varied. 
The lines are marked as for the /3-functions. 



quarks mass in IR units (left column) , the variation of the /3- function (middle column) , and 
the variation of the gamma function (right column) as we scan over a range of Tq which 
includes the blue and red dashed curves of Fig. ^ 

For X = 2 (top row) the quark mass (left column) is negative for a wide range of Tq 
between the zero mass solution (dashed red vertical line) and the critical value of Tq where 
the UV behavior of the solution changes such that the mass is no longer defined (solid blue 
vertical line). The mass curve oscillates as we approach the vertical blue line of critical Tq 
as depicted in Fig. but the oscillations lie too close to line to be resolved at the used 
scale of Tq. As x is increased to 3.9 (middle row), the zero mass solution is driven very 
close to the critical Tq. For x = 4.25 (bottom row), the solutions with negative quark mass, 
and in particular the one with zero mass, have disappeared. 
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In the the middle column we plot as the thin black curves the /3-function as Tq is varied 
over the range of the left hand plots with a constant step size. The blue thick curve is the 
limiting /3-function that always terminates at an UV fixed point and corresponds to the 
blue curve of Fig. For x = 2 all /3-functions with positive quark mass (black thin solid 
curves) are running, including the zero mass solution (red dashed thick curve), whereas all 
walking solutions have negative quark mass (thin dashed curves)^". The running curves are 
the lowest ones and all of them basically overlap. The dotted thin lines are the /3-functions 
with bounce-back behavior in the UV, which occur for values of Tq that are smaller than 
the critical one. Increasing x to 3.9 (middle row), the zero mass solution moves into the 
walking region, and approaches the critical one marked with thick blue curve. As x is 
further increased (bottom row), the zero mass solution disappears by joining the blue 
curve. 

Similar, but slightly more complicated behavior is seen for the gamma functions (right 
column). The dashing and coloring have the same meaning as for the /3-functions. Notice 
that 7/r asymptotes to —3 as A —)• for the zero mass solution (red thick dashed curve) 
whereas for the solutions with a finite mass it asymptotes to —1, as expected. The solutions 
with negative quark mass (thin dashed curves) have a zero of the tachyon, which shows up 
as a pole in 7/r. The main change in the plots as x increases from 2 (top) to 4.25 (bottom) 
is the movement of the solution with the UV fixed point (blue thick curve) towards smaller 
A, as all solutions it "crosses" change drastically. Otherwise the gamma functions (i.e., 
those solutions left of the blue curve) are roughly independent of x. 



Finally, we comment on the m — t- limit and discuss also Fig. 11 in this same limit. 



For X < Xc the background converges smoothly to the m = one with chiral symmetry 



breaking in this limit, as also seen from the /3- and 7-functions on Fig. 16 (top and middle 
rows). Taking m — )• for x < Xc in Fig. |TI| defines the enveloping curve of the fixed mass 
curves, which diverges for x = Xc- (We shall discuss the m = case in more detail in 
Sec. 10 below.) For x > Xc, the scale ratio Auv/^IR diverges for m — )• (see Fig. ^), and 
we can consider two different limits. If we keep Auv fixed as m — )• 0, the IR scale is driven 
to infinity and the background converges pointwise to the one with identically vanishing 



tachyon and an IR fixed point, discussed in Sec. g.5.1. If we instead keep the IR scale fixed. 



the background converges towards the one having an UV fixed point of Sec. D.2| (see the /3 



and 7-functions shown with the thick blue line on the bottom row of Fig. 16). Notice that 
the backgrounds with vanishing tachyon and x < x^ are not connected to the backgrounds 
having a finite quark mass by any limiting procedure, which is in line with our expectation 
that they are unphysical. Indeed we shall show in Section ^ that these solutions have larger 
free energy than the ones with nontrivial tachyon and chiral symmetry breaking. 

G. Extracting UV coefficients from numerical solutions 
G.l Extracting free energy differences 

As pointed out in Section ^, we need to check numerically which one of the two solutions 

^"As mentioned above, there are also positive quark mass solutions which basically overlap with the thick 
blue curve. Our resolution is not enough to resolve these. 
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with vanishing quark mass, the one with chiral symmetry breaking or without, minimizes 
the free energy. In order to do this we need to extract A (for fixed A) defined in Eq. (| 
for both solutions and then calculate the energy difference through Eq. ( |9.10 ). 



The corrections involving A are highly suppressed 



in the region which is under 



perturbative control (log(rA) is small). Extracting A directly from the numerical solutions 
of A and A is practically impossible, since it is difficult to require the two solutions to 
have the same A to a high enough precision. Therefore we study variations in X, which is 
invariant in scalings of r. In principle we could match the numerically extracted variation 
of X to the correction term in Eq. ( p. 14 ). This is doable (except for values of x very close 
to the critical one Xc — 3.9959), but a large uncertainty in the value of A still remains. 
Therefore we proceed as follows. We substitute Xq + Xi in the equation ( 8.19 ), where Xi 



is treated as a small perturbation. From the linearized equation, we can solve Xi exactly: 



Xi{X) = Xcexp 



l6V{X)Xo{X)yX + 3V'{X) + 3Xo{XfV'{X) 
6F(A)Xo(A)2 



dX 



(G.l) 



where Xc is a constant. So, if Xq is known, this equation gives Xi in Eq. ( 9.12| ) to all 
orders in log(rA). Since the UV expansion of Xq is also known, we can use that to expand 
Xi and to calculate the relation between the constants A and Xq by using Eq. ( |9.13 ). A 
straightforward calculation gives 



A 



WiA^X, 




log(9yi/8)(23yi2 -64^2) 



(G.2) 



16V{X)Xo{Xy/X + 3V'{X) + 3Xo{Xyv'{X) 



14 64^2 \ 
\^9yiA2^9A Wfx) 



dX + 



6V{x)XQixy 

32 1 



9Vi 



To obtain the free energy difference between two solutions, we calculate numerically the 
variation of X and match with Eq. ( p.lj ) , where we use either of the two solutions to evalu- 
ate Xq and the integral numerically. Since Xi given by Eq. ( p.lj ) is a good approximation 
already at small r (log(rA) does not need to be small) we can obtain the difference AXc 
between the solutions to a good precision. The value is then used to calculate numerically 

and AS through equations (|9.1C ) and ( p. 21) . 



G.2 Extracting the chiral condensate at m = 

Let us then discuss how the value of the chiral condensate can be extracted from a nu- 
merical solution to the differential equations. Recall that our method for constructing the 
backgrounds requires tuning the normalization of the tachyon in the IR (Tq) such that the 
quark mass vanishes. This procedure is however limited by the numerical precision of the 
solution, and therefore exactly zero quark mass cannot be obtained. As it turns out, this 
makes direct extraction of the condensate value difficult, since the approximately linear 
term ~ mr of the tachyon dominates over the cubic solution ~ ar^ in the deep UV. This 
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is particularly problematic as x approaches Xc, since the ratio of the UV and IR energy 
scales grows, pushing the asymptotic UV region to smaller r while the condensate value 
oc a decreases. 

Therefore, we employ a subtraction procedure. In the UV, the tachyon solution (with 
the mass rrii as small as can be achieved) can be written as 

ri(r) = miTmir) + aMr) . (G.3) 

The UV expansions of the two (approximately) linearly independent solutions T^^i — 
£r{— log(Ar))'-" and Ta-{r) ~ £r^{— log(Ar))~'-' are given in Appendix p. We will remove 
the remaining mass term by subtracting another solution 

T2(r) = m2Tm{r) + (T2r„(r) (G.4) 

which is calculated at the same value of x but for different Tq such that m2 is about the 
same order as mi. Since the mass terms dominate at very small r, it is easy to extract 
1712 /mi to a high precision, and construct'^^ 

Tsir) = Ti(r) - ^T2{r) = L - T^r) . (G.5) 

m2 V fn2 J 

If a is analytic at m = 0, 

ai = aW+a«m, + a(2)m2 + ... , (G.6) 



we obtain 



Tsir) 



^(0) [ 1 _ ) _^ o-(2)mi(mi - m2) + 



T.(r) . (G.7) 



If the masses rrii are sufficiently small, the quadratic correction oc a^'^^ as well as higher order 
terms can be neglected, and the difference Ts is proportional to a^^^ to a good precision. 
The reliability of the subtraction procedure can be tested by varying, say, m2 and checking 
that this does not affect the result. 

Now fj^'^) can extracted from Tg in a straightforward manner. However the result still 
has quite large error bars in particular for x close to Xc, because of numerical uncertainty in 
the deep UV. We do an additional trick to remove this problem. The extracted Ts extends 
close enough to the UV singularity to ensure that the tachyon is decoupled from A and A. 
On the other hand, we can obtain a decoupled tachyon solution with zero quark mass by 
solving first the background functions A and A with T = 0, inserting these in the tachyon 
EoM, and solving that by shooting from the UV. We can basically extend this solution as 
close to the UV singularity as we want. Matching the two solutions in the region where 
both are reliable, we can extend Ts easily up to ~ j--^oo ^j^g \JY. 



There is a subtlety here: the UV scales A which appear in Eq. (D.9) will be different for the two 
solutions. This needs to be fixed, e.g., by scaling one of the solutions such that the corresponding solutions 
for A match as A — >■ 0. 
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Finally, to extract we consider the expansion of the tachyon at small A: 



logr5(A)-log( 1-^) -log^~log(T(°)+logT,(A)-log£ (G. 
, <t(°) 8 /4/ii 39V?-64V2\, WiX 



which is obtained by using Eqs. ( D.9 ) and ( p.2| ). The coefficients Vi and hi were defined 
in Eqs. ( p.4D , and A = Auv is the scale of the UV expansions. We match this expansion 
with the extended Ts, and extrapolate to A = to obtain log(cj*^'''' /^^)- 

H. Details on BKT scaling 

In this Appendix we discuss some technical details on the BKT scaling. First, there is a 
subtlety with matching the approximation ( |10.11| ) with the tachyon solution in the IR (for 
which we were unable to write analytic approximations) . Since the approximation contains 
a periodic function, the boundary conditions at r ~ will not fix the tachyon uniquely: 
one can let the tachyon oscillate first, and do the matching only after the oscillations have 
created n zeroes, where n = 0, 1, 2, . . .. Therefore, in Eq. ( 10.15| ) we should actually have 
(including the solutions where the tachyon changes sign) 



Vk(A* - A,) log ^^7rn + 0{1) . (H.l) 
r 



An analogous result was found in [^] in a case that was analytically more tractable. As 
in their case, we expect that the solution with no nodes (n = 0) is the one that minimizes 
the free energy, as we already verified numerically in Section ^. 

One can actually do the matching procedure even more precisely, and in particular 
derive results for K and K. We need to study the approximation ( |10.11| ) more closely: 
actually it can be matched with the other solutions only in the vicinity of the zeroes of the 
sine function. At r = rjR we have 

nnR)=^Qpf— ) sin(V^(A,-A,)log^ + 0) (H.2) 



T'irm) ^ Cfp^ 2 sin ( V'^(A* - Ae) log ^ + 

+\/ k{K - Ac) cos (^\/k(A* - Ac) log ^ + cf)^ 

It appears that the cosine factor in the derivative can be neglected due to the small factor 
a/ k{\* — Ac) multiplying it. However, in this case the ratio T{riji)/T' {rm) = l/2riR is 
fixed so we cannot match both the function and its derivative with the (nontrivial) IR part 
of the tachyon solution. Therefore, we indeed need to be close to a zero of the sine. 



riR 



sm 



y^(A*-Ac)log ^ + i)= 0{^/x7^c) , (H.3) 
which brings in the required dependence of the solution on both Cfp and </>. 
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Matching the solution ( 10.11 ) towards the UV is only possible close to its node, as 
well^^. This gives 

4> = 0(\/A* - Ac) . (H.4) 
Consequently, the argument of the sine in the approximation ( 10.11| ) must change by a 



positive multiple of tt as we move from the UV to the IR, so that Eq. ( [H.lj ) is actually 
written as 

v/k(A, - Ae) log ^ ~ 7r(n + 1) . (H.5) 
r 

where n is the number of the zeroes of the tachyon solution^^. Notice that Eqs. ( |10.13| ) 
and ( |10.14| ) also receive extra factors from the refined scaling argument. We have now 



C{p ~ err 



UV 



\/A* — Ac 



njv 



1 



\/A* — Ac 



(H.6) 



instead of Eq. ( |10.14| ). The additional square root factors cancel in the result for a, and 
they could in any case be neglected as a subleading correction to the exponential scaling. 

In order to write down the results for K and K we recall that the definitions of Xc{x) 
and Xc read 

G(Ac(x),x) = 4 ; G{X,{x,),Xc)=4 , (H.7) 

respectively, where A*(x) was defined by Vg{X^:{x)) — xVj(^{X^:{x)) = 0. Near the critical 
point we have the two expansions 



d 



4 + — G(Ac,x)(A* - Ac) 



G(A=K, X 

where k = -^G{Xc, x). By using Eq. (H.5) we finally see that 



4+^G'(A,(x),x)|^.^^Jx-Xc) 



(H.^ 



r-uv 



where 



K 



TT 



exp 



TT 



K{n + 1) 
\/A* — Ac 



exp 



K{n + 1) 

V^c - X 



K 



TT 



AG{X,,x) J-iG{Kix),x)\ 



(H.9) 



(H.IO) 



and setting n = gives the result for the solution which has no tachyon zeroes and lowest 
free energy. 



^■^Proving this turns out to be more tricky than in the IR, since the solutions ( 10.10 ) and ( 10.11 ) have 
the same power behavior, and therefore they apparently join smoothly for large values of the sine function. 
However, this is a fake effect due to the roughness of the approximation ( lO.lC ): we have seen that in the 
UV the power law of the tachyon actually changes from 2 to 3 within a range of log r which is small with 
respect to l/y/X* — Ac. Therefore the same argument as in the IR applies, and the sine function needs to 
be 0{\/ A* — Ac) in order the matching to work. 

'^''Our arguments do not exclude zeroes of the tachyon solution near the endpoints of the validity of 
the approximation (10.11), i.e., near r — r and r = rm. However it is reasonable to expect that the 
approximation joins the UV and IR solutions in a smooth manner, so that no extra nodes appear, which 
can also be verified numerically. 
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Finally, let us notice the free energy scaling result for the solutions having several 



tachyon zeroes. Following the arguments of Sec. (10.2), we see that the energy difference 
between the solution with vanishing quark mass and the solution with n tachyon zeroes 
scales as 



AS 



4K(n + l) 
"v/A* — Ac 



UyCXp 



X 



(H.ll) 



Since we have numerically checked that AS is positive, this result verifies that the n 
solution has the lowest free energy for x — )• Xc. 
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